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Abstract 

We introduce various exchange interactions in a model Hamiltonian for 
Ce 3+ ions in cubic symmetry with three configurations (Z ,/ 1 ,/ 2 )- With the 
impurity pseudo spin Si = 1/2, our Hamiltonian includes: (i) One-channel 
S c = 1/2 Anderson model; (ii) Two-channel S c = 1/2 Anderson model; (iii) 
An unforseen one-channel S c = 3/2 Anderson model with a non-trivial fixed 
point; (iv) Mixing exchange interaction between the Tqj and the T$ con- 
duction electron partial wave states; (v) Multiple conduction electron par- 
tial wave states. Using the third-order scaling (perturbative renormalization 
group) analysis, we study stability of various fixed points relevant to various 
exchange interactions for Ce 3+ ions in cubic symmetry. 

PACS Nos. 74.70.Vy, 74.65. +n, 74.70.Tx 
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I. INTRODUCTION. 



Since the introduction of the orbitally nondegenerate Anderson model []IJ, realistic ex- 
tensions of this model have been considered @,[|. Hirst |2| described the Anderson model 
including the spin-orbit interaction and the crystal electric field effect based on the group 
theory. A detailed study of possible low temperature physics along this line was given by 
Nozieres and Blandin || for the magnetic ions in normal metal hosts (primarily transition 
metal ions). Nozieres and Blandin discovered the multichannel Kondo models which have 
non-trivial fixed points in the overscreened cases, but concluded that the realization of this 
model has a difficulty in real materials due to the channel asymmetry. Cox introduced the 
two-channel quadrupolar Anderson model |§] to explain the physics of UBei3, and listed 
possible crystal symmetries || for the existence of the two-channel exchange interaction in 
rare earth and actinide compounds or alloys. However, a detailed study is required for the 
stability of the two-channel fixed point against possible symmetry-breaking fields. Though 
in a different context, the metallic two- level system (electron-assisted atom tunneling in a 
double well) was shown to map into the two-channel Kondo model ||. 

In this paper, we study in detail a one-impurity Anderson model for Ce 3+ with three 
configurations in cubic symmetry in the same spirit as Nozieres and Blandin did for the 
transition metal ions. Our study can be generalized to any / atomic configurations in cubic 
symmetry. As we shall show below, the inclusion of the detailed energy structure in the 
atomic orbital leads to rich physics. In summary, we find various exchange interactions with 
the impurity pseudo spin Si = 1/2: (a) One-channel S c = 1/2 exchange interactions for 
r 6 7 conduction electrons (magnetic "Kramers'" doublet orbitals); (b) Two-channel S c = 
1/2 exchange interaction for Tg conduction electrons (quartet of orbitals); (c) One-channel 
S c = 3/2 exchange interaction for r§ conduction electrons; (d) Mixing exchange interactions 
between Tqj and Tg conduction electrons. In addition, (e) the multiple conduction electron 
partial wave states enhance the exchange couplings. In connection to the normal states of 
the cuprate superconductors |7j and the recent discoveries of non- Fermi liquid systems || in 
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heavy fermion materials, the relevant question is whether the ground states of the models 
(b) or (c) can be realized in the presence of other exchange interactions. 

We will answer this question using perturbative renormalization group (third-order scal- 
ing) arguments. We find three kinds of stable fixed points with the impurity spin Si = 1/2: 
one-channel Fermi liquid fixed points; two-channel non- Fermi liquid fixed point; and three- 
channel non-Fermi liquid fixed points. In addition, a "zoo" of unstable fixed points are 
discovered. At these unstable fixed points, all the possible exchange interactions can exist 
in the octuplet manifolds of the conduction electrons. Furthermore the multiple conduction 
electron partial waves generate the multi-channel unstable fixed points. Some of our fixed 
points possess dimensionless couplings which are rigorously in the perturbative regime, while 
for the above mentioned one-, two-, and three-channel cases this is not so. Nonetheless, since 
we are only after the existence of the fixed points, the scaling theory has proven qualita- 
tively accurate for non-trivial two- and three-channel fixed points, and the one-channel fixed 
points are well known to flow to strong coupling. Thus we have confidence in the correct- 
ness of our qualitative conclusions regarding stability of the various fixed points, provided 
we understand the modification for the one-channel case. 

Similar results and conclusions have been reached by Zarand || in the two-level system 
Kondo model (electron assisted tunneling of an atom in a double well) where he employed 
a large number of channels (large real spin) and thus rigorously controlled the perturbative 
renormalization group theory. We have verified the essential correctness of our scaling results 
for the S c = 3/2 conduction electrons with a numerical renormalization group study of the 
model (also combined with some conformal field theory analysis of the model) JT^JTTj] A 
recent conformal field theory study of single channel, large conduction spin Kondo models, 



with special emphasis on the S c = 3/2 case; the results are compatible with ours | I2]| . Since 
this fixed point is ultimately unstable against partial wave mixing exchange terms, it may 
be of little physical relevance. However, an interesting theoretical question is raised by our 
study of the fixed point, bolstered by an independent numerical renormalization group and 
conformal field theory study ||11|| . Namely, the universality class of the fixed point is not 



that of the three-channel Si = 1/2 Kondo model as conjectured by Tsvelik and Wiegman 
fl"3| and utilized in their Bethe-Ansatz treatment of the multichannel model, a point also 
noted by Zarand in his study of the novel unstable fixed points of the TLS Kondo model 



]. These ideas will be discussed more extensively in Ref. [Ill 



Our paper is organized as follows: Section II contains an overview of the generalized 
Anderson model Hamiltonian and the underlying group theoretical motivations for under- 
standing the various couplings which arise. Section III contains a precise and systematic 
discussion of the hybridization matrix elements which feed into the model together with 
derivation of the effective Schrieffer- Wolff couplings between conduction and 4/ electrons 
upon integration out of virtual charge fluctuations. Section IV contains the scaling analysis 
of various models which may be obtained in limiting cases of the full Hamiltonian including 
all possible conduction electron symmetry states without adding in multiple partial waves. 
Section V discusses the effect of including multiple conduction electron partial waves in 
each symmetry channel. Finally, we summarize and offer conclusions in Section VI. Three 
appendices detail finer points about the derivation of mixing matrix elements (App. A), the 
multiplicative renormalization group applied here (App. B), and the general form of the 
third order scaling equations for our models (App. C). 



II. MODEL HAMILTONIAN I: OVERVIEW OF ENERGETICS AND 

SYMMETRY. 

The model Hamiltonian we consider throughout this paper is the generalized Anderson 
model . Schematically, this Hamiltonian has three parts, defined by 

H = H cond + Hf + H hyb (1) 

where H con d defines a broad flat conduction band of width D, Hf defines the quasi-atomic 4/ 
states on the Ce site as modified by metallic screening of the direct Coulomb interaction and 
the crystalline electric field (CEF), and H hyb = H i+Hi 2 describes the hybridization between 
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the 4/ and conduction electrons. Hqi describes mixing of f° — f 1 excitations with conduction 
electrons, and H u describes mixing of J 1 — / 2 excitations with conduction electrons. We shall 
describe Hh y b in detail in the next section; in this section we shall focus on the symmetry and 
energetic aspects of the first two terms. In particular, we shall emphasize the symmetries of 
tensor operators from the conduction and / sectors which may be coupled together in the 
effective exchange interactions mediated by virtual interconfiguration excitations. 

Before proceeding further, we provide a brief physical description of the different irrep 
labels of the cubic group: 

(1) T 6 and T 7 are magnetic Kramers' doublet (Kramers'). Two degenerate states can be 
transformed into each other with time reversal symmetry operator on them and will be 
labeled with pseudo-spin j or j. That is, the T 6 and T 7 CEF states are similiar to the 
J = 1/2 angular momentum manifold. 

(2) T 8 is a magnetic quartet (r 8 = T 3 ® T 7 ). With the time reversal symmetry operator, 
the four- fold degeneracy can be decomposed into two pairs. These two pairs are disjoint in 
terms of time reversal and will be labeled by orbital indices, ±, which correspond to the 
quadrupolar or orbital T 3 degree of freedom - stretched (3z 2 — r 2 like) or squashed (x 2 — y 2 
like) local orbitals guaranteed degenerate by the cubic symmetry. The T 7 degree of freedom 
will be labeled as j and | as above. Overall, the r 8 quartet is labeled by ± j / |. Under 
certain circumstances it is more favorable to view the T 8 manifold as an effective S — 3/2 
manifold, as we shall discuss below. 

(3) The orbital singlets r\ and T 2 do not need any further label. The T 2 is a pseudoscalar, 
the lowest even parity version being (x 2 — y 2 )(y 2 — z 2 )(z 2 — x 2 ). 

(4) The orbital (non-magnetic, or non-Kramers') doublet T 3 will be labeled by ± as in the 
above. In this case, time reversal symmetry operator does not transform one state to the 
other as for a Kramers' doublet. 

(5) The magnetic triplets T 4 and T 5 are labeled by 0, ±1. Under the time reversal symmetry 
operation, the ±1 states transform into the =Fl states, respectively. On the other hand, the 
state transforms into itself. That is, the T 4 and T 5 CEF states are similiar to the J = 1 
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multiplets under time reversal symmetry operation. 

In the absence of energetic considerations, there are 14 (91) different degenerate eigen- 
states in the f 1 (f 2 ) configuration. These states can be progressively split by the strong 
on-site Coulomb "Hund's rule" interactions (/ 2 ), spin-orbit coupling, and the crystalline 
electric field (CEF). The energy hierarchy in Cerium is, roughly, that the interconfiguration 
splitting are several eV, the exchange splittings are about 0.2 — 0.5 eV for the f 2 , the spin 
orbit splittings are order 0.3 eV, and the typical CEF splittings are from a few to tens of 
meV. Specifically, the / — f 1 interconfiguration excitation energy is E{f 1 ) — E(f°) ~ — 2eV, 
and the f l -f 2 energy difference is E(f 2 )-E(f v ) « 4eV 0. Traditionally, E(/ 1 )-E(/ ) is 
referred to as tf while E(f 2 ) — E{f 1 ) is referred to as 6f + U, meaning the effective Coulomb 
interaction is about 6 eV. 

For the f 1 configuration, the energy level structure is simple due to the absence of the 
on-site Coulomb interaction. Spin-orbit coupling first results in J = 5/2,7/2 multiplets. 
In most Ce 3+ compounds and alloys, this energy splitting is always about 0.3 eV |L5| . 



Furthermore the crystalline electric field will split each J = 5/2,7/2 multiplet. The CEF 
can mix two different J's but with the second order perturbation correction of order | < J = 
5/2\Hcef\J — 7/2 > \ 2 1 A LS which is of the order of 0.1 x A CEF . Neglecting the mixing of 
two different J's in the presence of the CEF, the irreducible representations (irreps) of the 



full rotation group, D 5 / 2 ,D 7 / 2 decompose in the cubic field as |L6| 



d 5/2 = r 7 © r 8 , (2) 
d 7/2 = r 6 © r 7 © r 8 . (3) 

Here and T 7 are the magnetic Kramers doublet's irreps and Tg is the magnetic quartet 
irrep of the cubic double group. The eigenstates of r 7 8 for J = 5/2 multiplets are given in 



Table I [ ITU . We will remove the superscript J = 5/2 below when the context is clear. The 
eigenstates of r 6i7;8 for J = 7/2 multiplets are listed in Table |. The choice of the overall 
phase is arbitrary in defining the CEF eigenstates. The above convention will make the 
effective exchange interaction a simple form. Throughout this work we shall assume that 



the J = 5/2, T 7 doublet lies lower than the J = 5/2, T 8 quartet. 

For the f 2 configuration, the existence of the strong on-site Coulomb interaction compli- 
cates the energy level structure. The LS Russel-Saunders scheme applies to the / electron 
in Ce ||18|| . First in the absence of both spin-orbit coupling and CEF, total orbital angular 
momentum L and total spin angular momentum S are good quantum numbers, and thus 
energy eigenstates can be sorted out by these quantum numbers. In the presence of spin- 
orbit coupling, the total angular momentum J = L + S is a good quantum number. There 
are seven possible values of J starting from through 6. When Ce 3+ ions are sitting in a 
cubic environment, each J multiplet further splits into CEF energy levels (see Table |J). For 
example, the f 2 Hund's rule ground multiplet J = 4( 3 if 4 ) decomposes as 

At = Ti e r 3 © r 4 e r 5 . (4) 

All in all, there are 7 IVs, 3 ^'s, 9 T^s, 9 r 4 's, and 12 T^s in the f 2 configuration. Here 
Ti and T 2 are singlet irreps, T 3 the non-Kramers doublet irrep (nonmagnetic), and r 45 
the magnetic triplet irreps of the cubic point group. The eigenstates of r 1)3j4 5 for J = 4 
multiplets are in Table |TTT| 

We now consider the conduction electrons. According to the Anderson model picture, the 
conduction electrons can hop on and off the atomic orbitals at the impurity site. The Bloch 
electrons may be projected at the impurity site into three irreps in the cubic symmetry: 
T 6 , Ty, r 8 . For the / shell in Ce, we expect the I = 3 conduction electron partial waves 
are most strongly coupled to the / shell. For the isotropic hybridization, only the I = 3 
components can hybridize with the / orbitals. In the presence of the spin-orbit coupling, 
the total angular momentum is a good quantum number. Thus we have J = 5/2,7/2 
conduction electron partial waves, these J multiplets further split into CEF irreps in crystal 
environments as mentioned above. When the real charge fluctuations are removed from the 
model system in the Kondo limit, we have to construct the tensor operators representing each 
CEF state for the f 1 configuration and the projected conduction electron CEF states. The 
relevant tensor operators are for re,7,8 CEF states. We list the decomposition of conduction 
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electron operators into irreducible tensor opertors in Table |V| using the character table [16 



In the direct product, the first CEF states are written as ket, and the second as bra. The 
1^6,7 tensor operator (2x2 tensor) is a direct sum of a charge operator (Ti) and a pseudo spin 
operator (T4). Indeed, the Schrieffer- Wolff transformation leads to two interaction terms: 
the spin exchange interaction and the pure potential scattering term. The relevant T 8 tensor 
operators are the 21^ triplets for our model. In the conduction electron Tg tensor space, 
one of the two IYs gives rise to the pseudo spin S c = 1/2 operators with two degenerate 
orbital channels and the other to the pseudo spin S c = 3/2 operator with one channel. The 
interpretation of the orbital channels is that these are shape degrees of freedom of the local 
T 8 states-one is "stretched" or 3z 2 — r 2 like, and the other "squashed" or x 2 — y 2 like. Thus, 
the r 6; 7 conduction electrons provide two different channels with the pseudo spin S c = 1/2, 
while the T 8 conduction electrons provide two degenerate channels with S c = 1/2 or one 
channel with S c = 3/2. 

As a further complication, mixing is possible between r 6 7 and T 8 conduction electrons. 
The mixing tensor operators can be read off from Table [TV]. That is, each mixing can provide 
one pseudo spin representation. Hence from the group theory, we conclude that there are 
in total 6 possible exchange interactions between the f l J = 5/2 T 7 pseudo-spin and the 
three-symmetry conduction electrons. We are going to derive these in the next section. 



III. MODEL HAMILTON! AN II: ANDERSON HYBRIDIZATION. 

In this section we derive explicitly the forms of the hybridization matrix elements between 
excited f°, f 2 states and f 1 states which are mediated by conduction electrons hopping on 
and off the impurity. We also derive the contributions to the effective exchange interactions 
obtained by integrating out virtual charge fluctuations to the f°, f 2 states. 

The conventional theory of Ce 3+ impurities assumes an infinite on-site Coulomb interac- 
tion, which removes the f 2 configuration from consideration. When we relax our assumption 
about the strong on-site Coulomb interaction and we include the detailed atomic energy 
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structure in cubic crystals, we can construct a variety of model Hamiltonians. In this work 
we assume that the magnetic doublet f l J = 5/2T 7 states lies lower than the f 1 J = 5/2T 8 
quartet states. The one-channel S c = Sj = 1/2 Anderson hybridization is present between 
f° and f x J = 5/2 T 7 states. This gives rise to the conventional single channel Kondo 
coupling. From the hybridization between f 1 and f 2 configurations, various exchange inter- 
actions arise from virtual charge fluctuations, which include: 

(i) the one-channel exchange interaction for magnetic doublet Tqj conduction electron or- 
bitals; 

(ii) the two-channel S c = 1/2 exchange interaction for Tg conduction quartets; 

(iii) the one-channel S c = 3/2 exchange interactions for T 8 conduction electrons; 

(iv) mixing exchange interactions between the T 8 and the r 6 7 conduction electrons; 

(v) multiple conduction electron partial wave states. 

Below we list all the possible terms arising from the Anderson hybridization in the 
cubic symmetry. All the possible selection rules for the hybridization are listed in Table 
|V[ Since the Anderson hybridization is irreducible tensor of r 1; the selection rule is solely 
determined by the three components: the conduction electron partial wave states and the 
other two neighboring atomic configuration states. In the f 1 configuration, we only consider 
the ground CEF states f 1 .] = 5/2 T 7 . In the f 2 configuration, there are all 7I\, 3r 2 , 9r 3 , 
9r 4 , and 12r 5 CEF states. 

For definiteness, we will consider isotropic hybridization. By this we mean that the 
original hybridization potential is taken as isotropic and this will be projected down to the 
appropriate irreps of the cubic point group symmetry at the Ce ion. Then J c = 5/2,7/2 
conduction electron partial waves are mediating the tunneling process between neighboring 
configurations of the atomic orbitals. 

The mixing between the f° and the f 1 configurations is simple and straightforward to 
evaluate, since the f° is a singlet with Ti character. This is given by the Hamiltonian terms 

= ^EE 4m c \f° >< f 1 ; 3 C m c \ + h.c. (5) 

e j c m c 
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in the absence of the CEF, or 

Hoi = V J2 E 4r cdc \f° >< f'-JcTM + h.c. (6) 

in the presence of the CEF. Here the conduction electron Bloch states are symmetry adapted 
into appropriate partial wave states labeled by cubic irrep indices r c at the impurity site. 
These projected quantum labels play an important role in determining the ground state 
properties of the model system. 

However, the mixing between the f 1 and the f 2 configurations becomes complicated due 
to the fractional parentage and the Clebsch-Gordan coefficients. First we consider the case 
in the absence of the CEF. 

H?2 ] = V EE E E Hjcm^fhrmlfLSJM) 

e jcm c j\m\ LSJM 

xc\ Jcmc Ifjimi >< fLSJM\ + h.c. (7) 

Here we use the reduced matrix elements are defined 

HfjimrJcmclfLSJM) = <j 1 m 1 | / ic77lc \f 2 LSJM > 

= K( 3c ; fMfLSJ) < 3c m c ;j imi \JM > . (8) 

In the last line we used the Wigner-Eckart theorem. The prefactor K(j c ; f 1 ji\f 2 LSJ) is 
the fractional parentage coefficient (see the Appendix |A| for more details). Though we 
projected the atomic Fock space using the symmetry quantum labels, these states are not 
energy eigenstates of the atomic Hamiltonian. The Coulomb (exchange) interaction will 
mix multiplets with the same J. This will introduce another unitary transformation from 
the symmetry states \f 2 LSJM > to the energy eigenstates. These energy eigenstates still 
respect the total angular momentum and are labeled by JM with extra energy quantum 
label. We eschew this complication in favor of working in the simpler LS coupling limit, 
which is accurate for Ce ions in any case. 

As will be discussed below, the cubic Clebsch-Gordan coefficients for different J's are 
equal to within a prefactor for f 1 T 7 CEF states. Hence the mixing between the different 
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multiplets with the same value of J will not change the matrix form (hybridization or mixing 
matrix) . 

In passing we note that the same approach was undertaken for the f 2 and f 3 configu- 
rations including only the Hund's rule ground multiplets ||19|| . In principle, the discussion 
we give here for the f 1 , J = 5/2, T 7 doublet should equally well apply to the f 3 , J = 9/2, T 6 
doublet considered in Ref. [19|], which however contends that only a fermi liquid state is 
possible. We shall defer a discussion about this different conclusion between the two works 
to the last section of the paper. 

We now turn to a systematic discussion of the hybridization matrix elements between f 1 
and f 2 and the effective exchange interactions mediated by virtual f 2 fluctuations. 



A. The hybridization between f l J = 5/2 T 7 and / 2 ri,r 2 . 

As shown in Table |V|, only the r 7 (r@) conduction electrons can mix with the f 1 T 7 to reach 
the f 2 Ti(T 2 ) states. The corresponding cubic Clebsch-Gordan coefficients can be written in 
a very simple form 

< Jc I>; ^T 7 (3\J 2 T 1 > = {-lf 2 - a (p,a) < j c T 7 ; >, (9) 

< j c T Q a-h 7 (3\J 2 T 2 > = (-l)Va-a < j c Y 7 -h 7 \J 2 V 2 > . (10) 

The phase dependence derives from time reversal symmetry. The common factors, which are 
independent of the CEF degeneracy labels, < j c r 7 ; |r 7 | J 2 Ti > and < j c r 7 ; |r 7 | J 2 T 2 > 
are listed in Table |VT| and |V11| . 

Then the Anderson hybridization can be written as 

ffiaCTx) = VY, E ("ir 1/2 K(j c ; A \ f 2 LS : J) < Jc T 7 ; |r 7 | JT t > 

j c LSJ 1 1 

x4 jc r 7a If^ra >< f 2 LS : JTi| + h.c, (11) 
H 12 (T 2 ) = VY, E (-l) a - 1/2 K(j c ; f\ | fLS : J) < j c F 6 ; |r 7 |jr 2 > 

ea j c LSJ A A 

xir 6a l/^IVa >< fLS : JT 2 \ + h.c. (12) 
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In the mixing, the degeneracy index for a time reveral pair is active, i.e., can change. In 



the local moment limit, the Schrieffer- Wolff transformation j20] leads to the one-channel 
antiferromagnetic exchange interaction 



Hi(Ti) — ^2 Jij S c r r ij(0) ■ Sr 7 , Jij — yJiJj, (13) 

ij=5/2J/2 



Ji= jc = ^ j E2{ LSJ)- El lK{3c] fl l 1 fLS : J)|2 1 < JcT7] l T7lJTl > |2 ' 



(14) 



and 



ffi(r 2 ) = J5 c r 6 (0).Sr 7 , (15) 

J = V , \K(-; f- | fLS : J)| 2 | < ~r 6 ; -T 7 \JT 2 > | 2 . 

f£jE 2 {LSJ) - Ei K 2 ,J 2 lJ n 1 2 '2 1 1 

(16) 

Here the conduction pseudo-spin is S c = 1/2. We also introduced the energy levels for atomic 
states in the f 1 and the f configurations. E\ is the energy level for the f 1 J = 5/2 T 7 state 
and E 2 (LSJ) for the f atomic state with quantum labels of orbital (L), spin (S), and total 
angular momenta (J). 

B. The hybridization between fj = 5/2 T 7 and / 2 r 3 . 

Only the Tg conduction electrons can mix between the f l J = 5/2T 7 and the / 2 r 3 's. The 
relevant cubic Clebsch-Gordan coefficients are 

<i c r 8 nV;^r 7 a|J 2 r 3 n> = (-l) a+1 / 2 (a', «) (n» < j c T 8 ; ^T 7 \ J 2 T 3 > . 

(17) 

The overall phase factor derives from time reversal symmetry. The common factor 
< jc^; |r 7 | J 2 r 3 > is listed in Table |V11J|. Hence the hybridization can be written as 



H 12 (T 3 ) = V £ £ J2 (-l)- 1 / 2 K(j c ; f\ | fLS : J) < 3c T 8 ; |r 7 | Jr 3 > 

xir 8 „„ If-IV* >< / 2 LS : JI>| + /i.e. (18) 
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In the hybridization, the time reversal pair index a =|, \ is active while the orbital index 
n — ± does not change and is not active. Thus the time reversal pair (f / j) gives rise 
to the pseudo-spin and the orbital doublet (±) provides two channels. The Schrieffer- Wolff 
transformation leads to the two-channel exchange interaction 



Hi 



Jij S c r 8 ijn(S)) ' Sp 7 , Jij 



J, 



l=Jc 



ij=5/2,7/2 n=± 

LSJ E 2 (LSJ) - Ei 



(19) 



\K(j c ; f\ | f 2 LS : J)| 2 | < j c r 8 ; ^T 7 \ JT 3 > | 2 . 



(20) 



Here the conduction pseudo-spin is S c = 1/2. 



C. The hybridization between f 1 J = 5/2 T 7 and / 2 r 4 . 

The T 7 or T 8 conduction electrons can mix between the f 1 J = 5/2T 7 and the / 2 r 4 's. 
The hybridization has the following form 

#i 2 (r 4 ) = ^EE Kti* f \ I / 2 £S = J) < JcTs; |r 7 | jr 4 > 



J2a V2 



,t 

ej c r 8 -a 



l/ 1 ^ >< f 2 LSJT$\ 



x < 



x < 



(21) 



+ ( 4r 8 -T I/ 1 t> ~VS 4 r8+i I/ 1 I> ) < f 2 LSJT 4 l\ 

^ +( 4rw If 1 i> -Vs 4 r8+T I/ 1 T> ) < / 2 ^jr 4 l| 
E K{j c -f\\fLS-.j) < Jc r 7 ;|r 7 |jr 4 > 

E«ir ra l/ 1 ^ >< / 2 LSJr 4 0| 

+V2 4 rrT I/ 1 1>< / 2 LSjr 4 i| 

^ +^4^1 I/ 1 i>< PLSJY A l\ _ 
The hybridization looks more complicated compared to those for r 1; T 2 and T 3 CEF states. 
However, this complication comes from our special labeling of the CEF degeneracy as 
we shall show below. We list the numerical common factors < jj^s] |T 7 | J 2 r 4 > and 
< j c r 8 ; |r 7 | J 2 T 4 > in Table [TX|. Note that the explicit dependence on the degeneracy 
labels is left out in extracting the prefactors. 
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The Schrieffer- Wolff transformation leads to a new exchange interaction model with the 
S c = 3/2 conduction electrons interacting with the impurity pseudo spin Si — 1/2 in the T 8 
conduction electron space 



Hi — JijS c r 8 ij(0) ■ Sr 7 — J <S> S'cr^O) • Sry, — y JiJj, 
2\V\ 2 
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\K( 3c ; f'l I f 2 LS : J)| 2 I < Jc r 8 ; \v 7 \JY A > | 2 . 



(22) 



(23) 



Here S c = 3/2. The exchange coupling can be rewritten in the matrix form. The correspond- 
ing representations of the angular momentum operators are (|+ |>, |+ |>, |— |>, |— |>) 



St 



s: 





/-3 





\ 


1 





3 





2 





1 







V 





-1/ 


/ 





\ 





V3 











-2 


Vv^ 





/ 


(° 








Vs\ 

















V3 








U 





-2 


/ 



-3 > 



2° 



S* 



v 7 ^ 
-2 



5* 



-v 7 ^ 
-V3 -2 



(24) 



We can show that these matrices transform into the canonical SU(2) form after the rear- 
rangement of the states: 

,3, 1. 



(1+ T>, 1+ l>, I- T>, I- !>) - (^(-^ >> 1 2 2 > '~ ][ \\ >; 



;)>)• 



This extra complication comes from our biased labeling of the CEF degeneracy toward the 
two-channel exchange interaction. /,From a viewpoint of the two-channel exchange interac- 
tion, the induced exchange interaction can be written as 
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Hi — J <S> (— r° - 2t 1 ) <S> 5*p g S* 7 , 



7- = — — r — 



Vs 1 



1 



(25) 

a/3 

2 I"*, r« = --r 3 + — r 1 , r 2 = r 3 . (26) 

Thus we can see that the exchange interaction generated by the f 2 T 4 CEF states are harmful 
to the two-channel Kondo effect. 

On the other hand, the T 7 conduction electrons are coupled to the impurity spin with 
the ferromagnetic exchange coupling 

Hi — — JijS c r 7 ij(0) ■ Sr 7 , Jij = \J Ji-Jj, (27) 

rl5 



J, 



= E 



(28) 



Here = 1/2. In this case, the exchange interaction, mediated by the f 2 T 4 states, de- 
grade the one-channel exchange coupling which is ever present between the /° and the f l 
configurations. 

Another complication arises through the / 2 r 4 CEF states. Two conduction electrons 
of F 7 and Tg can mix with each other. Since it can be easily generalized to the multiple 
conduction partial waves, we present the model for the one partial conduction wave case. In 
the T 7 © r 8 conduction electron space, all the possible exchange interactions are 



Hi = J 1 

j = J 7 k° + J 8 r° + J 1s k 1 + J 88 r\ 



(29) 
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x < 



D. The hybridization between J 1 J = 5/2 T 7 and / 2 r 5 . 

The T 6 or T s conduction electrons can mix between the f 1 J = 5/2T 7 and the / 2 r 5 's. 
The hybridization has the following form 

h 12 (t 5 ) = V£ E K( ]c -f l \\fLS:J) < Jc r 8 ;|r 7 |jr 5 > 

E a v^cj,^^ \fa >< f 2 LSJT 5 0\ 

+ ( V3 4r 8 _ T I/ 1 T> +4r 8+ l I/ 1 l> ) < / 2 ^/r 5 l| 
^ +( y/Z 4^ I/ 1 |> +4 r8+T I/ 1 T> ) < f 2 LSJT 5 l\ 

K(j c] f l \ | fLS : j) < Jc r 6; |r 7 |jr 5 > 

-Ea4 r6a l/ 1 ^ >< / 2 L5Jr 5 0| 
WHcT.t I/ 1 1>< / 2 mr 5 i| 

l/U></ 2 ^Jr 5 lj ^ 

The Schrieffer- Wolff transformation leads to a new exchange interaction with the S c = 3/2 
conduction electrons interacting with the impurity pseudo spin Si = 1/2 in the T 8 conduction 
electron space 



x < 



> . 



(32) 



Hi — ^ JijS c r s ij(fy ' 'Stt, Jij — \J JiJj, 

ij 

2\V\ 2 



J - = E TO*; /I I ^ = ^)l 2 1 < ^s; ^r 7 | > | 2 . 



(33) 



(34) 



Here S c = 3/2. This conduction electron pseudo spin is different from that defined in the 
previous section. The corresponding representations of the angular momentum operators 
are (|+ |>, |+ |>, |- T>, |- |>) 
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o -Vs o o 
V o o o o j 

We can show that these matrices transform into the canonical form after the rearrangement 
of the states: 

(1+ !>,!+ I>,|- t>,|- I» - C-l|| >,l§(~) >,-!§<-§) >,l§§ ». 

/From a viewpoint of the two-channel exchange, the induced exchange interaction can be 
written as 



E x = J <g> (-r° + 2r l ) <g> S» g S* 7 , 



1 ^ 1 3 , i 2 

r — _ T T — — —j- -\ — _ T T = r 

2 2 2 2 



(36) 
(37) 



Thus in the Tg conduction electron space, there are three competitors: two-channel Kondo 
effect and the other two S c = 3/2 Kondo effect (overscreened) . 

On the other hand, the T 6 conduction electrons are coupled to the impurity pseudo spin 
with the ferromagnetic exchange coupling 



Hi — — ^2 JijScT 6 ij{fy ' Sr 7 , Jij — y JiJj, 



2|Vf ] K(j c ; | fLS : J)| 2 | < Jc T 6 ; h 7 \JY A > ' 2 



LSJ 



E 2 (LSJ)-E 1 



(38) 



(39) 
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Here S c — 1/2. In this case, the exchange interaction, mediated by the / 2 r 4 states, degrade 
the one-channel exchange coupling. 

Another complication arises through the f 2 T 5 CEF states. Two conduction electrons of 
T 6 and T 8 can mix with each other. In the r 6 £g>r 8 conduction electron space, all the possible 
exchange interactions are 



Hi = f 



7 



O'i oi 
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(40) 
(41) 



(42) 



(43) 



E. The Tq © 1^ © Tg conduction electron space. 

In cubic symmetry, when we keep the lowest lying /'J = 5/2 T 7 CEF states, the most 
general exchange interaction form to each conduction electron partial wave species is 



h 1 = j*® % si,, 

J* = -h 7° + -h^ + J 8 r° + J 6 s Y + ^78 «* + Jss r\ 



(44) 
(45) 
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(46) 



(47) 



(48) 



The matrix algebra for 7, r is listed below. We note that in the subspaces, T 7 © T 8 or 
r 6 ©r 8 , the exchange interactions are self-contained. No new interactions are generated with 
renormalization. Also in the whole r 6 © T 7 © T 8 conduction electron space, the exchange 
interactions are self-contained. 



7V 



7 +2 r +2T, 



(49) 
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7* + 7 J = ~l k (cyclic), 

7V + 7V = l\ 

7V + 7 i 7 i = -7 - i r° + r k (cyclic) 

7V + r*7° = 0, for i = 0, x, y, z, 

r 7 + 7 r = 7 , 

r V + 7V = 7 fc (cyclic) , 

7V' + r j Y = 7 fc (cyclic), 

7 r + r 7 = 7 . 

r r — t , 

t t + r r = zr , 

r l + r j = -r k (cyclic), 

^xv -u^z zx 1 ^0 2 

T T — r y r = r r = r r , 

2 2 

^y^z _z-y xz 1 1 2 
r a r = t t = t t = — r H r . 

2 2 

+ /t- 7 = -K k (cyclic), 

K { K j + = -k° - \r°-T k (cyclic) 
k°t 1 + t 1 k° = 0, for i — 0,x, y, z, 

r l K j + /cV = -K k (cyclic), 

K l T j + T j K % = -K k (cyclic), 

reV + tV = -k\ 
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reV + 7^° = o, 
7 V + ^7° = 0, 



x/3 



7 V + ^7* = -yX, 
«V + 7*«* = 0, 



7°X + X7° = X- 



(72) 
(73) 
(74) 

(75) 
(76) 
(77) 
(78) 



X 



(79) 



Here x is given by 

/0 1 0\ 
10 

\fl 0/ 

Though this new term (%) arises in the intermediate calculation step, this term cancels out 
in the final result. 

In deriving the scaling equations, the following relations are useful 



2^8 ^88 2 ^ 68 ^ 78 ^ 



J 1 J' + J j r = (2J 6 2 - J 6 2 8 ) 7° + (2J 2 - J 2 8 ) k° + 

+J 68 ( 2 J 8 s — J& — Js ) l k + J78 ( — 2 Js8 — J7 — Js) K k 

+ ( ^68 — ^78 ~~ 2J 8 J88 ) T^- 



(80) 



Here (i,j,k) is cyclic, and 



J 1 J 1 = ( Jl + J 6 2 8 ) 7° + ( J 2 7 + 4) ^ + 



^8 + ^88 + o ("^68 + "^i 



78; 



+ ^68 ( ^88 + J6 + Js ) Y + J78 ( ~ J88 + J 7 + J 8 )K i 
1 



+ 



2 J?,Ji 



88 



(J 6g J 7g ) 



(81) 
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F. Summary 



In summary, we discovered various exchange interactions. In particular, a new type of 
exchange interaction is discovered: the S c = 3/2 conduction electrons interact with the 
impurity pseudo spin Si = 1/2. All the possible exchange interactions in our model study 
are summarized in Fig. [I]. 

The T 6 conduction electrons (pseudo spin S c = 1/2) are coupled to the impurity pseudo 
spin Si = 1/2 with the antiferromagnetic exchange coupling through f 2 T2 states, and with 
the ferromagnetic exchange coupling through f 2 T 5 states. The mixing between T 6 and T 8 
conduction electrons is also present. That is, the Tq conduction electrons can change their 
symmetry states (Tg) after scattered off the impurity pseudo spin. 

The T 7 conduction electrons (pseudo spin S c = 1/2) are coupled to the impurity pseudo 
spin Si = 1/2 with the antiferromagnetic exchange coupling through f 2 Ti states, and with 
the ferromagnetic exchange coupling through / 2 r 4 states. The mixing between T 7 and r 8 
conduction electrons is also present. That is, the I\ conduction electrons can change their 
symmetry states (r§) after scattered off the impurity spin. 

The Tg conduction electrons can generate two different conduction electron pseudo-spins: 
S c = 1/2 with two degenerate orbital channels and and S c = 3/2 with one channel. The 
T 8 conduction electrons are coupled to the impurity pseudo spin in a two-channel exchange 
interaction through / 2 r3, and in a one-channel exchange interaction through f 2 T^ states. 
Both models lead to the overcompensation of the impurity pseudo spin Si = 1/2. The 
two-channel case has been studied extensively using several techniques and leads to a non- 
Fermi liquid fixed point. Using the numerical renormalization group and the conformal field 
theory, the one-channel S c = 3/2 case is also shown to lead to the non- Fermi liquid fixed 



point [11 



We can give some plausible argument why some of exchange couplings are ferromagnetic 
while the others are antiferromagnetic. This argument follows from the addition of two 
angular momenta. As noted before, the CEF manifolds in the f 2 configuration can be 
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roughly evisaged as angular momentum multiplets. 

To hop into the orbital singlet r\ (T 2 ) in the f 2 configuration, the r 7 (r 6 ) conduction 
electrons have to be in the opposite pseudo-spin state compared to the / 1 T 7 (impurity) 
pseudo spin state. Thus the coupling should be antiferromagnetic. 

Now we consider the orbital doublet T 3 in the f 2 configuration. This orbital doublet is 
also represented as J = singlet in terms of time reversal. In this case, the relevant r 8 
conduction electrons are represented by pseudo-spin S c = 1/2 and should be in the opposite 
pseudo spin state compared to the f l Y 7 (impurity) pseudo spin state. Thus the coupling 
should be antiferromagnetic. ^From this argument, we can clearly see why the T 8 conduction 
electron cannot be represented by pseudo-spin S c = 3/2 for the virtual charge fluctuation 
to the T 3 CEF states. 

When we consider the magetic triplets T 4 and T 5 in the f 2 configuration, the exchange 
coupling can be ferromagnetic or antiferromagnetic depending on which CEF symmetry 
electrons hop into the atomic orbitals. These magnetic triplets can be represented by the 
J = 1 manifolds. For the T 6 or the T 7 conduction electrons which are represented by pseudo- 
spin S c = 1/2, the coupling should be ferromagnetic to form the J = 1 triplet state from 
the addition of S c — 1/2 and Sj = 1/2. As we noted before, the r 8 conduction electrons 
can be represented in both ways: S c = 3/2 with one orbital channel or S c = 1/2 with two 
orbitally degenerate channels. Hence depending on the spin representation, the coupling 
has different sign. With S c = 1/2 representation, the coupling should be ferromagnetic 
as explained above. This picture in fact agrees with our results that the virtual charge 
fluctuations to the triplets T 4 and T 5 in the f 2 are harmful to the two-channel Kondo effect 
for the r 8 conduction electrons. With S c = 3/2 representation, the coupling should be 
antiferromagnetic which can be deduced from the angular momentum addition of S c — 3/2 
and Si — 1/2 which has to be J = 1. 
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IV. SCALING ANALYSIS I: VARIOUS FIXED POINTS 



In this section we analyze our various model exchange interaction terms using the third 
order scaling argument - perturbative renormalization group (RG) 0.21-23]. At tempera- 
ture T, only the conduction electrons (thermally excited) inside the band of order T with 
respect to the Fermi level play an important role in determining physical properties. Thus 
we can integrate out the band edge states (virtually excited states) to find the effective 
Hamiltonian. Though the following analysis is restricted to the perturbative regime (weak 
coupling limit), we can derive qualitatively correct results out of this. For quantitative 
results, a full numerical renormalization group (NRG) study is required. 

Our strategy is as follows: we start with each exchange interaction term and identify the 
relevant fixed points. After that, we include other interaction terms step by step to see if 
the relevant fixed points remain stable. Finally, in the next section, we consider the effect 
of the multiple conduction electron partial wave states on the various Kondo effects. 

In summary, we find that three types of fixed points are stable in the presence of the full 
exchange interactions: 

(i) one-channel Fermi liquid fixed points (Tq and Tj conduction electrons); 

(ii) two-channel non- Fermi liquid fixed point (r 8 conduction electrons); 

(iii) three-channel non-Fermi liquid fixed points (some linear combination of r 8 and T 6 or T 7 
conduction electrons). In addition to these stable fixed points, we find a "zoo" of unstable 
fixed points at which various exchange interactions we can imagine are realized. When 
we consider the multiple conduction electron partial waves, the relevant stable fixed points 
remain stable with enhanced bare exchange couplings and new unstable multichannel fixed 
points are generated. 
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A. Structure of the fixed points 



Before discussing the scaling equations for various exchange interactions, we wish to 
establish some mathematical structure of the fixed points generated by the third order 
scaling equations. We consider the following generic exchange interactions 

H 1 = J l ® Si(0) Sj. (82) 

Here S c = Sj = 1/2 are spin operators and J 1 is the exchange coupling matrices whose 
dimension is 4 x 4 in our model case. We first consider some arbitrary SU(2) angular 
momentum operators 

U = M l ® - a\ (83) 
2 

Here o % are Pauli matrices and M % are undetermined at this point. Now we demand that 
L % satisfy the SU (2) spin algebra, that is, 

[ L\ V } = ie ijk L k . (84) 

Here is Levi-Civita antisymmetry tensor. Then we find the condition for M l 's to satisfy 

{ M\ M j } = 2 M k . (85) 

Here (i, j, k) is cyclic. Up to this point, all is a pure mathematical fun. Now we are going 
to infuse some physics into the above algebra. 

In the scaling analysis (the perturbative renormalization group), we remove the high 
energy states progressively (renormalize) to find the effective Hamiltonian which describes 
the same low temperature physics of the original Hamiltonian. When the effective Hamilto- 
nian does not change further with renormalization (removal of the high energy states), it is 
called the fixed point Hamiltonian. For the above general exchange interactions, the scaling 
equations up to the third order diagrams of Fig. || are (see the Appendices & 0) 

dg k 



dx 



g l g ] + g J g l \ -~ A g k ^ [g l g l + g'g J \ ■ (86) 
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Here k) is cyclic. All our exchange interactions we are going to study satisfy the cyclic 
property such that Tr[ ] is independent of the x,y,z indices. Then we can see that the 
equations for the fixed points have the same structure as in the 577(2) spin algebra. That 
is, the exchange coupling matrices have the SU(2) spin structure at the fixed points 

g* = g* M\ (87) 

g* = -4,-, (88) 
TrM 2 

A = ^ M 2 . (89) 

The fixed point can be further specified by A which measures the conduction electron spin 
size. The number of channels is determined by the ratio of the degeneracy of the relevant 
manifold and the spin size. Below we are going to illustrate this point explicitly. 

In passing we note that the same approach was independently applied to the metallic 
two- level system ||. 

B. The one-channel S c = 3/2, Si = 1/2 exchange interaction 

As shown in chapter II, the r 8 conduction electrons are coupled to the impurity pseudo 
spin Si = 1/2 with their pseudospin S c = 3/2 through the virtual fluctuations to the Z 2 ^^ 
triplet states. 

For comparison, it is well understood that the one-channel S c = Si = 1/2 exchange 
interactions (Tq and conduction electrons) lead to the Fermi liquid strong coupling fixed 
point, while the two-channel S c = Si = 1/2 exchange interactions (r 8 conduction electrons) 
give rise to the non- Fermi liquid non-trivial fixed point. 

The third order scaling equation is 

% = g 2 - \ s c (s c + l)(2s c + 1) g\ g = N(0)J. (90) 
ox 3 

Here s c = 3/2. We can find the stable fixed point g* = 1/5. 

We now argue that this new exchange interaction can lead to a non-Fermi liquid fixed 
point. In Fig. |3], we schematically show how the screening of the impurity spin can happen 
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by the conduction electron spin cloud. With the antiferromagnetic coupling, the energy is 
minimized when two conduction electrons of S* = —3/2, —1/2 are coupled to the impurity 
spin Sj = 1/2. This picture does not violate the Pauli exclusion principle. Thus the effective 
impurity spin is "SV" = 3/2. In the next shell screening, only the = 3/2, 1/2 can come 
closer to the effective impurity spin due to the Pauli exclusion principle. Hence the effective 
exchange interaction is antiferromagnetic. From this argument, we can see that the strong 
coupling fixed point cannot remain stable. Since the weak coupling fixed point is unstable, 
as we shall show, an intermediate fixed point must exist |§ . This is a non- Fermi liquid fixed 
point. We are separately studying this simple new exchange interaction using the numerical 
renormalization group(NRG) and conformal field theory flTT| . 

We are next going to show that this new exchange interaction competes with the two- 
channel Kondo effect. 



C. The Tg conduction electron space 

In the r 8 conduction electron space, two different spins can be realized as shown in 
section III. In summary, three different exchange interactions are possible via the f 2 T^^^ 

Hi{Ts) = h E Scn(O) • S F7 = J 3 r° <g> 5 C (0) • S F7 , (91) 

n=± 

HxiTt) = J 4 L c (0) -S r7 = J 4 ( - r° - 2r l ) <g> Sj(0) S* T , (92) 
Hi(T s ) = J 5 L' c (0) ■ S Tr = J 5 ( - r° + 2r J ) ® Sj(0) S* rr (93) 

Here Sr r is the impurity pseudo spin for the f l J = 5/2 Tf state. L c (0) and L c (0) are 
two different conduction electron pseudo spin representations for the same S c = 3/2 state 
manifold. Here we have rewritten the S c = 3/2 operators from a view point of the two- 
channel exchange interaction. 

All combined in the r 8 conduction electron space, we find 

H 1 = J i ^> Si S\ (94) 
j = J 8 r° + J 88 t\ (95) 
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Js — J3 — Ja — J5, J88 — J5 — Ji- (96) 

Physically, the r° term leads to the ordinary two-channel exchange interaction while the r % 
terms are channel-symmetry breaking and channel-mixing. Then we can derive the scaling 
equation for the matrix g 1 = N(0)J l in a compact way. In the r 8 conduction electron space, 
the third order scaling equations are 

^ = 9* ~ \ 9 2 s 8 ~ 9s [ 9 2 8 + g 2 88 I 98 = N(0)J 8 > 0, (97) 

(><M - 288 [ gt + gls ], 988 = iV(o) J 88 . (98) 



dx 

We can identify the two-channel fixed point (<7g, #g 8 ) = (1,0) and two other interesting fixed 
points (3$,^,) = (-1/5, ±2/5) 

g* = I M, (99) 
M* = -r° ± 2 t\ (100) 
A = §•§,». (101) 

corresponding to the exchange interaction models generated by / 2 r4 i 5 in this restricted 
space, respectively. The linear analysis shows that all these fixed points are stable. Solving 
the above scaling equations, we can construct the flow diagram (see Fig. |j). There are three 
linear separatices along which all the exchange interactions derived purely from the / 2 r 3 , T 4 , 
or T 5 CEF states lie. Along the g$ axis, the two-channel S c — Sj — 1/2 exchange interaction 
due to virtual f 2 T 3 charge fluctuations is located. Along the line g 88 = 2g 8 , the one-channel 
S c = 3/2, Si = 1/2 exchange interaction due to virtual f 2 T^ charge fluctuations is located. 
Along the line g 88 = —2g s , the one-channel S c = 3/2, Si = 1/2 exchange interaction due 
to virtual f 2 T 3 charge fluctuations is located. ^From this flow diagram, we can get the 
following qualitative results: (1) When 2g$ > |g 88 |, two-channel fixed point ground state is 
realized; (2) When g 88 > and 2g 8 < |g 88 |, the new fixed point ground state (generated 
by / 2 r 5 ) is realized. (3) When g 8S < and 2g 8 < |g 88 |, the new fixed point ground state 
(generated by / 2 T 4 ) is realized. 
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From a viewpoint of the two-channel exchange interaction, the above model Hamiltonian 
has channel asymmetry and exchange coupling anisotropy without the r 1 terms. From a 
well-known scaling argument 0,^5|,^1J , the model will flow to the one-channel fixed point (the 
channel with larger z component of exchange coupling) in the absence of the channel-mixing 
term. Though the r 3 term breaks channel symmetry, the additional term from r 1 mixes the 
two different channels. That is, the channel asymmetry due to r 3 term is apparently washed 
out by the channel mixing term restoring to the two-channel fixed point in the two-channel 
Kondo effect parameter regime. 

We note that the S c = 3/2 fixed points have a net coupling strength (along the separat- 
ices) of g* = 1/5, safely in the perturbative regime so that the scaling analysis is trustworthy. 



D. One-channel & two-channel exchange interaction 



When we keep only the three atomic states f°, f 1 J = 5/2 T 7 , and / 2 r 3 , we find the one- 
channel (r 7 ) & two-channel (r 8 ±) Anderson hybridization interactions with the impurity 
pseudo spin Si = 1/2. This simple model shows competition between the Fermi liquid fixed 
point of the one-channel Kondo effect and the non- Fermi liquid fixed point of the two-channel 
Kondo effect [[24[] . For completeness, we will include the scaling analysis of this model here. 
We studied this simplified model using the non-crossing approximation . 

It can be deduced from the scaling theory that the low temperature and the low energy 
physics is dominated by the one-channel or the two-channel Kondo effects depending on 
their relative magnitude of the antiferromagnetic couplings. With the introduction of the 
exchange coupling matrix, it is more convenient for the derivation of the scaling equations 
to rewrite the one-channel and two-channel exchange interactions in the following form 



H x = J ® S 




Si, 









\ 


J = 





■h 












Js) 



(102) 
(103) 
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— * — * 

Here S c and Si are S = 1/2 operators. The scaling equations of our simple model Hamilto- 
nian up to the third order diagrams are 

|| = g 2 -igTr[g 2 ]..., (104) 
g = iV(0)J. (105) 

The scaling equations in components are 

|| = 9r - I 97 [ 97 + H }, 97 = N(0)J 7 > 0, (106) 
|| = 9l ~ I 9s [ 9 2 7 + 2g 2 8 ), g 8 = N(0)J 8 > 0. (107) 

Here x = log(D/T). We can identify three fixed points related to one-, two-, and three- 
channel Kondo effects. The one-channel, strong coupling fixed point (g 7 ,gs) — (oo,0) is 
beyond the perturbative regime, but is stable leading to the singlet ground state and Fermi 
liquid spectrum p7|. The three-channel fixed point (2/3, 2/3) is stable along the line g-j = g$ 



in the gj — g$ plane, but unstable for any small perturbation from gj = g$. Finally, the two- 
channel fixed point (0, 1) is stable leading to the logarithmically divergent thermodynamics 
at zero temperature. ^From the scaling analysis, we can infer the ground state physics: 
one-channel for Jj > J§; two-channel for Jj < J§; and three-channel for Jj = Jg. 



E. Complete © Tg exchange coupling space 

In the J c = 5/2 T 7 © r 8 conduction electron subspace allowing for other than virtual 
/ 2 T 3 charge fluctuations, all the possible exchange interactions can be written in a simple, 
compact form 

H x = f ®S£(0) Si r , (108) 
J 1 = J 7 k° + J 8 r° + J 78 k' + J 88 t\ (109) 

The matrices are defined in the section III. We already studied this model in some specific 
cases in the previous sections. We can find the scaling equations up to the third order (see 
the Appendices & 0) as 
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% = 9r - \ 978 ~\97[g 2 7 + 2(g 2 78 + & + <&) ], (no) 



<9#78 



+ ^(07 + 0s) ) 



X - g 78 [ g 2 + 2{g 2 78 + g 2 8 + g 2 88 ) ], (111) 



^ = gl - \<ik - \<$% -\g%[<$ + % 7 2 8 + gl + gh) ], ( 112 ) 

= -#8#88 - ^re - ^ #88 [ 97 + 2(^78 + #8 + #8s) ], ( 113 ) 

g 7 = N(0)J 7 , g 78 = N(0)J 78 , g 8 = N(0)J 8 , g 88 = N(0)J 88 . (114) 

From these scaling equations, a linearized analysis reveals that some of the fixed points we 
found in the previous sections remain stable. The relevant fixed points (g 7 , g 78 , gl, g 88 ) are: 

(i) (oo,0,0,0): one-channel fermi liquid fixed point; 

(ii) (0,0,1,0): two-channel non-fermi liquid fixed point remains stable against any small 
perturbation away from this fixed point 

g* = M, ikP = r , (115) 
A = i.5 T o (116 ) 

As discussed previously, A measures the conduction pseudo-spin size and the number of 
channels. Here the conduction electrons are described with the pseudo-spin S c = 1/2 and 
with two orbitally degenerate channels (Trr = 2); 

(iii) (2/3,0,2/3,0): three-channel non-fermi liquid fixed point 

g* = I M, (117) 
jkP = /t° + T , (118) 
A = ^U° + r°). (119) 

is stable against the — r 8 mixing interaction or the T 8 ± channel mixing interaction, but 
unstable for any perturbation away from g 7 — g 8 ; 

(iv) (0, 0, -1/5, ±2/5): S c = 3/2 fixed points are 

gl = I Ml, (120) 
M l ± = -t° ± 2 t\ (121) 
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A = - • - r°. (122) 
2 2 K J 

The g+ fixed point is stable, but the other S c = 3/2 fixed point g*_ becomes unstable 
against the T 7 and T 8 mixing exchange interaction (J78). As we shall show below, the 
mixing interaction between T 6 and T 8 conduction electrons is relevant to the fixed point. 

There are another two unstable fixed points which we did not consider before. The fixed 
points are g* = 2/11, g* 8 = -2/11, g* 78 = 0, g* 88 = ±4/11 

M i = k° - r° ± 2 t\ (124) 

A = \ ■ Ik + 1 4 r°. (125) 
2 2 2 2 v ; 

At the fixed points the reduced exchange interaction is 

H x = l. (§M + S' c (0)) -Sr 7 . (126) 

Here S c — 1/2 and = 3/2. That is, two different conduction "bands" with different spins 
are coupled to the impurity spin. 

All the fixed points we considered so far have g? s = 0. Since the I^ — T 8 mixing interaction 
is relevant perturbation for some fixed point, we might imagine that some fixed points with 
a finite value of g 78 are possible. There are in total 6 fixed points of which only two are 
stable 

2 2 4\/2 4 

^0? = -9> 9l = - 9 , <4 = ±^p & = -g, (127) 

B:g; = l> 9l = ~\, & = ±^, %*8 = -^ (128) 

C : £7 = - — , dl = — , ^78 = ±— , £ 88 = — — ■ (129) 

Out of these, only the type A fixed points are stable. The type B & C fixed points are 
unstable against a perturbation along one principal direction. 
The A fixed points are 
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3 T ' 



A 



(130) 

(131) 
(132) 



The fixed points can be represented with the direct sum of three S c = 1/2 conduction 
electron pseudo spins. We further note that the complex coupling matrices can be simplified 
using the orthogonal transformation in the conduction electron channel space 

( -\ \ /-l 0\ (\ \ 



M x = 



0-1 

V o -l o j 



M y = 



1 

V o 1 oj 



M z = 



0-10 
\0 -I) 



(133) 



After rearrangement and redefinition of the states with some phase and another orthogonal 
transformation to diagonalize the spin-flipping terms in the channel space, we can rewrite 
the exchange interaction at the fixed points in a concise form 

3 



(134) 



n=l 



This is the three-channel S c — Si — 1/2 exchange interaction. We further note that the 
following exchange interactions are all equivalent and result in the same physics, because 
the SU(2) algebra does not change (we call this the chiral - "handedness" - symmetry). 



H 1 = J S c (0) ■ Sj, 

H X = J{- S X C S^ - S»fi? + S z c Sf ), 
H 1 = J{ S*S* - S^Sf - S z c S z j ), 
H 1 = J( - S X C S* + S y c S y j - S z c S z j ). 



(135) 
(136) 
(137) 
(138) 



We can show the equivalence between them by redefining the states with the introduction 
of additional phase. 

Type B & C fixed points can be written in a standard form 
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Sb = \M, (139) 

M < = |«°-|^±^«*-ir*, (140) 
A= 1 -2 ( k° + t° ). (141) 

The spin representation is S c = 1. 

gc = ^ M, (142) 
AT = -| «° + | x°± ^ «* - | t*, (143) 
A=^(«° + r°). (144) 

The spin representation is S c = 5/2. 

In the r 6 ©r 8 conduction electron space, we get the same physics as in the above analysis. 
The difference is that the S c = 3/2 fixed point stable against the T 7 — T 8 mixing becomes 
unstable against the r 6 — T 7 conduction electron mixing. Hence both of the S c = 3/2 fixed 
points cannot remain stable in the presence of the mixing interactions between the r 8 and 
the r 6j7 conduction electrons. 

F. The Tq(BT 7 (B T 8 conduction electron space 

In this section, we are going to study the full exchange interactions possible in our model 
approach. The extension of the scaling analysis to the full interactions is straightforward. 

#i = J' <g> Sj(0) S l Fr , (145) 
X = J 6 7° + J 7 k° + J 8 r° + J 68 7 * + J 78 k< + J 88 t\ (146) 

Third order scaling equations are 

^ = 3l-\ Sk -\gs[gl + 97 + 2(<? 6 2 8 + g 2 78 + gl + g 2 88 ) ], (147) 
^ = 97 ~ \ 9i8 - I 97 [ gl + g 7 + 2(4 + g 2 78 + g 2 8 + g 2 88 ) ], (148) 

2 ~^ 2 ~^ r 2 2 i ~^ r2 2 /2 2 2 2\i 

■^7 = #8 - 2^88 - 1 [ 9 68 + 978 \ ~ g ^ 8 t #6 + #7 + %68 + #78 + #8 + #8 8 ) ], 

34 



%!8 

dx 

dx 

dgss 
dx 



068 



= ~978 



988 + g^ + tfs) ) 



2 ^68 [ gl + g 2 7 + 2((? 6 2 8 + (? 7 2 8 + gl + :1 



\ 978 [ gl + g 2 7 + 2{g%8 + 07 2 8 + 8 2 J ' 



.If 2 2 1 1 2 2 2 2 

+ 2 L 068 - 078 J - 2 88 i 06 + 07 + 2 l068 + 078 + 08 ' 



g 6 = N(0)J 6 , g 7 = N{0)J 7 , g 78 = N(0) J 78 , ^ = iV(O) J 8 , (? 88 = iV(O) J 8 



(149) 



(150) 



(151) 



(152) 



(153) 



We can find various fixed points of the above scaling equations. We already considered the 
cases (1) g 78 = g* 68 = 0; (2)g* 8 ^ 0,g* 8 = 0; (3) g* 78 = 0,g* 68 ^ 0. One relevant question is 
to ask if the two-channel fixed point of the T 8 conduction electron manifolds remains stable 
in the presence of the full interactions. The two-channel fixed point turns out to be stable. 
There arises the four-channel fixed point due to the expansion of the conduction electron 
space 



W = 7° + k° + r°, 
A = I.|( 7 ° + « + r ). 



(154) 
(155) 
(156) 



This fixed point as well as the three-channel fixed points remain stable against a small 
perturbation of the mixing interactions (J68,^7s), but they are unstable away from g 6 = 
g 7 = g 8 . There is also an unstable two-channel fixed point when we fine tune g% = g 7 



g* = M, 
M 



7° + k°, 



A = i-f(7° + «°). 



(157) 
(158) 
(159) 
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This fixed point remains stable against the conduction electron mixing interactions but 
becomes unstable away from g e = g 7 . 

The two S c = 3/2 fixed points are unstable in the presence of the full mixing interactions 
between the r 6j7 and the T 8 conduction electrons. Since the g^s and gjs terms are relevant 
perturbations to the S c = 3/2 fixed points, we may imagine that these fixed points flow to 
the other fixed points with a finite value of mixing exchange couplings. We shall address 
this question shortly. 

Another relevant question is to ask if the new type A three-channel fixed points discussed 
in the previous section remain stable. These two fixed points which have 

&a = \m, (160) 

M i = -\{K» + T»)± 2 ^,e- 2 -T\ (161) 

A=^U° + r°) (162) 

are stable, modulo the separate equality of exchange constants in the T 6 or T 7 space with 
those of the r 8 space. 

Due to the expansion of the conduction electron space (r 6 ), other fixed points related 
to the types A, B, C (in the Y-j © r 8 conduction electron space) considered in the previous 
section arise 

Sl = \M, (163) 

M* =7° + r° )±^ «* -\t\ (164) 

A = ^(7° + -° + r°). (165) 

= I M, (166) 
M* = 7 ° + ^ -^t°±^^-^, (167) 

A=^7° + l-2(.° + r ). (168) 
S*c = ^ M, (169) 
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i^ = 7 o_| B » + ^ T » ± ^«'_| T ', (170) 

All these fixed points are unstable. 

We now consider the case: g 78 ^ 0, (7g 8 7^ 0. There are 8 new fixed points which are not 
discussed in the previous sections 

g* = ^ M, (172) 
M^ = 3 1 °- 7 -K°- 7 -r°±2V3Y ±^5|c* + lr«, (173) 



M* = -- 7 + 3 K ° - \ r° ± ^ 7 * ± 2^3 «' - \ r\ (174) 
A= 7 -. 9 -( 1 ° + k + t°). (175) 

These 8 fixed points are all unstable. The spin representation is S c = 7/2. 

In summary, various fixed points have been identified. See Table [XI] for a list of all 
the possible fixed points found from the third order scaling equations. The stable fixed 
points are: (i) one-channel fixed points; (ii) two-channel fixed point; (iii) three-channel fixed 
point. We also identified a "zoo" of unstable fixed points. At these unstable fixed points, 
various exchange interactions are possible: (iv) the higher conduction electron pseudo spins 
(S c = 3/2,5/2,7/2) are interacting with the impurity pseudo spin Sj = 1/2; (v) the multi- 
channel conduction electrons with different size of pseudo spins interact with the impurity 
pseudo spin. The higher conduction electron pseudo spin models overscreen the impurity 
pseudo spin and thus can lead to non-trivial non-Fermi liquid fixed points. 



V. SCALING ANALYSIS II: MULTIPLE CONDUCTION PARTIAL WAVE 

STATES. 

Up to now, we have not included the multiple conduction electron partial wave states. 
By this we mean, e.g., additional T 7 states of odd spatial parity about the Ce 3+ site which 
may mix with the principle J = 5/2 T-j states. In this section we are going to study the 
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effect of them on various Kondo effects. The point is that the multiple partial waves lead to 
the enhancement in the exchange coupling of the original exchange interaction model. 

In the cases of simple one-channel or two-channel exchange interactions, the algebra is 
very simple. When various exchange interactions are studied at the same time, the third 
order scaling equations become complex, large matrix differential equations. As an example, 
we analyze the matrix scaling equations in the r 8 conduction electron space. 



A. Two-channel exchange interaction model. 

We consider the two-channel exchange interaction with some arbitrary number (N) of 
Tg conduction electron states. Then the Anderson hybridization interaction reads in the 
two- channel sector 

#1 = EEt" 1 )^ 172 ^ c lr lna l/ lr V* >< / 2r >l + h.c. (176) 

e ina 

We focus on one / 2 T 3 CEF state for the moment. Here the index i labels the different r 8 's of 
the conduction electron partial waves projected at the atomic site. Then the Schrieffer- Wolff 
transformation leads to the effective exchange interaction 

^l = EI! J iJ ScF s ijn(0) ■ Sr 7 , (177) 
ij n=± 

ScT s ijn{0) = X C Il>a 2^ C fT 3 8 n/3' ( 178 ) 
en' 

i 2V 2 

Jij = y/JiJj, Ji = — ■ (179) 

e 2 — ei 

Here ei(e 2 ) is the energy of the atomic state f 1 T 7 (f 2 T 3 ), respectively. Note that the orbital 
degeneracy index n in the T 8 manifold is not active in the exchange interaction and thus 
provides the two degenerate channels. Except for the unphysical case of the same form of 
the projected conduction electron DOS for each T 8 symmetry, we can not diagonalize both 
the projected conduction electron Hamiltonian and the above interaction at the same time. 

We demonstrate below how the physics can be simplified in this unrealistic, simple case: 
the same form of the projected conduction band DOS for each T 8 symmetry. J^- S c r 8 ij n 
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can be diagonalized in the T 8 label space with the orthogonal transformation. Since the 
conduction electron part is a unit matrix in the T 8 label space, the diagonalization in the T 8 
label space can be done at the same time for the conduction electron part and the exchange 
interaction part 

Hi = J E*SU°)-Sr 7 , (180) 

n=± 

Scn(0) =J2J2 C tna ^a/9<V n/ J, (181) 
ea e >/3 Z 

C «* = S C ^na, (182) 

i 

We note that only one eigenvalue is non-zero and is given by J. All the other eigenvalues are 
identically zero. We have not included the other transformations for the conduction elec- 
trons which are irrelevant. The other components of the transformed conduction electrons 
(nonbonding combinations) are decoupled from the impurity spin. We conclude that the 
multiple T 8 conduction electrons give rise to the enhanced two-channel exchange coupling 
with J = J2iJi an d the nonbonding combinations decouple. 

We now consider more general case with all the excited f 2 T 3 CEF states considered. In 



this case, the exchange coupling is symmetric = Jji but without the relation J^- = y JjJ,-. 

— * 

In general, the simplification as above does not occur when we diagonalize Jy S c r 8 ij n - After 
diagonalization, we may enumerate the eigenvalues of in the descending order from the 
maximum positive value (Ji) to the smallest value (Jn). That is, J\ > J 2 > • • • > Jn- Then 
we can write the interaction in the diagonal form 

N 

#i = EE Ji S* n (0) - S Tr . (184) 

i=l n=± 

If there are any negative J's, their magnitude will be renormalized to zero at zero temper- 
ature. If Jj's are nondegenerate, as generically expected, the model system will flow to the 
two-channel fixed point for the most strongly coupled set of conduction electron states. If 
J\ = J 2 7^ J3, the system will flow to the four-channel fixed point, and so forth. 
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Though we have already diagonalized the exchange couplings in the T 8 conduction elec- 
tron label space for the same projected DOS, we can give the scaling argument for more 
general forms of both the symmetric exchange couplings and the projected DOS 

#i = EE 4; SUn(O) • £, (185) 

ij n=± 

J ij = J Jl , J«>0. (186) 

The index n labels the orbital degeneracy in the T 8 manifold and the indices i,j label the 
conduction electron partial waves. In the scaling analysis, it is more convenient to introduce 



the exchange coupling matrix = ^/A^(0)A^(0) Jy. Here iVj(e) is the projected DOS for 
the i-th T 8 conduction electron. Then the third order scaling equations can be written in a 
very compact form (g is an N x N matrix) 

|| = g 2 - g Trg 2 . (187) 

Now we can diagonalize the exchange coupling matrix and we assume it is already done. 
Then we can write the scaling equations for each diagonal component Gi 

*| = G?- Gi £Gj. (188) 

Whenever the diagonal elements are positive, they will grow initially, e.g., with decreasing 
temperatures. However, the largest G will determine the low temperature physics. There 
are several fixed points of which only one fixed point is stable according to third order scaling 
analysis. Third order scaling analysis says that the fixed point {G\ — 1, G\ — G\ — • • • — 
G* N = 0) is stable. This fixed point is none other than the two-channel fixed point with the 
initial enhanced exchange couplings. 

In general, the 2M-channel fixed point is G\ = G* = ■ ■ ■ = G* M = 1/M, G* = for 
i > M + 1 or 

1 (Im 0\ 

* = s( o) « 189 » 

In the fixed point coupling matrix, Im is an M x M unit matrix. This 2M-channel fixed 
point has one eigenvalue of — 1/M(irrelevant), M — 1 eigenvalues of 1/M (relevant), and the 

40 



remaining eigenvalues of —1/M (irrelevant, corresponding to G* = for i > M + 1). This 
2M-channel fixed point is unstable along the M — 1 principal directions and stable along 
the other remaining principal directions. Hence this 2M-channel fixed point is stable only 
when the bare couplings satisfy the relation: G± = ■ ■ ■ = Gm > Gi for i > M + 1. In the 
M = 1 case, this fixed point becomes stable as mentioned above, since the unstable principal 
directions are absent. 

In the case of the one-channel exchange interaction models (I"V and Tq conduction elec- 
trons), we get the same results as in the two-channel exchange interaction model case. 



B. Multiple Tg conduction electron partial waves. 

Now we study the effect of multiple T 8 conduction partial waves on the Kondo effect in 
detail allowing for the S c = 3/2 fixed point 

#1= [j 8 ®T + J 8 8®T i ] ®Si(0)S\ (190) 

Then it is straightforward to find the scaling equations up to third order (see the Appendices 
i&0) 



dx - Ss 2 



^ g 8 ^ 8 -g 8 Tr{g 8 2 + g 8 y, (191) 
1 +g88g 8 ]-g88Tr{g2 + g2 8 }. (192) 



dx 2 

Here g 8 and g 88 are reduced exchange coupling matrices. In general, when we consider the 
multiple conduction electron partial waves, the exchange couplings become matrices. 

For definiteness, we first consider two T 8 conduction partial waves, say, coming from 
J c = 5/2, 7/2 multiplets. In this case, the exchange coupling matrices are 2x2 matrices 

(193) 

To solve the above scaling equations we may introduce the following decomposition of the 
exchange coupling matrices 
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/Jll 


Jl 2 \ 




(Jn 


Jn\ 






; Jss — 






W21 


J22) 




U21 


J22) 



gs = 9o P° + 9i P 1 + 9s P 3 , (194) 
gss = 9o P° + 9i P 1 + 93 P 3 - (195) 

Here p° is a 2 x 2 unit matrix, and p*'s are Pauli matrices. These matrices are defined in the 
space of the conduction electron partial waves where the T 8 manifolds are derived. Since the 
exchange coupling matrices are real and symmetric, there is no p 2 term. Then the scaling 



equations in components read 

d -§ = \9\ 2 -\\9\ 2 -29o{\9\ 2 + \9\ 2 ), (196) 

^ = 2g ogi - Mi - 2 9l ( \g\ 2 + \~g\ 2 ), (197) 

^ = 2g g 3 - g g 3 ~ 293 ( \g\ 2 + \g\ 2 ), (198) 

^ = - g .~g-2g (\g\*+\g\*), (199) 

^ = ~\( 9o9i + 9o9i )-2g 1 ( \g\ 2 + \g\ 2 ), (200) 

^ = ~\( 9093 + 9093 ) - 2g 3 ( \g\ 2 + \~g\ 2 ), (201) 

\9\ 2 = 9 2 o + 9 2 + 9 2 3, (202) 

\g\ 2 = 9 2 o + 9 2 i+9l (203) 

g ■ g = goh + gm + 9393- (204) 

There are several fixed points. First of all, the weak coupling fixed point g§ = gg 8 = is 
unstable. The four-channel fixed point 

gs = ^ ( V §88= f ) (205) 
2 \0 1/ V° 0/ 

is stable along ^o-axis, while unstable away from g\ = g 3 = 0. The two-channel fixed point 



gl = 7:(p + sm9 p' + cosO p 3 )^ [ . (206) 

2 v ' \0 0j 

Here 9 is an arbitrary real number. We get the diagonal matrix along the principal directions 
after an orthogonal transformation. This fixed point is stable with enhanced two-channel 
exchange coupling. 
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There is another class of the fixed points when g| 8 = ±2g|. The fixed point = 
-1/10,^=^ = or 



1 



1 0' 



8S = L . ( 207 ) 



10 



1 



is unstable, which corresponds to the two-channel S c = 3/2 exchange interaction. The 
relevant stable fixed points are g$ = — 1/5, \gl\ 2 + \gl\ 2 = 1/25 or 



1 



1 + cos sin \ 1 / 1 ~ 



§8 = "7 - "7 h (208) 

5 \ sinfl l-cos0/ 5 \0 y 

which are none other than new S c = 3/2 non- Fermi liquid fixed points with enhanced 
exchange coupling after the orthogonal transformation at the fixed points. 

We can generalize the above analysis with arbitrary number of conduction electron partial 
waves (N). We can determine the stability of the fixed points relatively easily. The case 
gg 8 = is equivalent to the case of the two-channel exchange interaction with multiple 
partial waves. There are several multi-channel fixed points of which only the two-channel 
fixed point is stable. In this case, the fixed points g§ are solutions of 

g 8 2 -Ag 8 = 0, A = Trg 8 2 . (209) 

The trivial weak coupling fixed point gg = is unstable. To find non-trivial fixed points, we 
introduce the orthogonal transformation to diagonalize the above matrix equation. Then 
we find the general 2M-channel fixed points (after the orthogonal transformation) 

1 (Im 0\ 

8s = M I - -I ' = ( 21 °) 







with A = 1/M . Here Im is an M x M unit matrix. M can have any value of M — 1, 2, • • • , N. 
Now we consider the stability of the fixed points. We will show that only the fixed point for 
M — 1 is stable. Close to the fixed points, we can linearize the scaling equations to find 

^ = g*5g 8 + <5g 8 g* - A 5g 8 - g* Tr{g*5g 8 + 5g 8 g*}, (211) 
^ = -\ [g 8 *5g 88 + <5g 88 g*] - A 5g 88 , (212) 
A = Tr{g*g*} = ^. (213) 
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Diagonalization of the first equation can be done by the same orthogonal transformation 
introduced for finding the fixed points. Here it is enough to consider a small perturbation 
along the principal directions, since any other forms of perturbation can be written as a linear 
combination of those along the principal directions. Then along the principal directions, we 
find 

/Mo* 1 2 M 

^ = _ Sgi - A; T SgL % = 1, 2, • • • , M, (214) 
dx M M 2 ~ ' ' ' ' K J 

^ = -j^ 5 9l i = M + l,M + 2,---,N. (215) 

Thus this 2M-channel fixed point is stable along the N — M principal directions (corre- 
sponding to % = M + 1, M + 2, • • • , N). And it is stable along one direction and unstable 
along the other M — 1 directions after the diagonalization for % = 1, 2, • • • , M. For a small 
perturbation along the principal directions of gg 8 = (in fact, those of gg), we find 

^T = -^88, i = l,V..,M, (216) 
-77**4, i = M + l,M + 2,---,N. (217) 



dx M 

Thus the 2M-channel fixed point is stable for any small perturbation of gss- We conclude 
that only the two-channel fixed point for M — 1 is stable when we consider the multiple 
conduction electron partial waves. 

Another class of fixed points gg g = ±2gg are analyzed in the same way. The identification 
of the fixed points are simple: there are M-channel S c — 3/2 fixed points of which only the 
case M — 1 is stable. The M-channel S c — 3/2 fixed point matrix is (after diagonalization) 

gss = ±2gs, (218) 
1 (Im 0\ 

^-5m( j- < 219 » 
Here Im is an M x M unit matrix. Now we study the stability of this fixed point using the 
linearized scaling equations close to the fixed point 
<9<5g8 1 

-g^T = gf^gS + SgS g8 - o [g88*g88 + ^g88 gs 8 ] - A *g8 
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-g* Tr{g*5g 8 + 5g 8 g* + g* 8 5g 88 + <5g 88 g* 8 }, (220) 

= ~\ [§88^88 + fegss + §8^88 + ^g 88 g 8 ] ~ A 5g 88 

-g* 8 Tr{g*<5g 8 + 5g 8 g* + g* 8 5g 88 + 5g 88 g* 8 }, (221) 
A = Tr{g*g* + g* 8 g* 8 } = ^, (222) 

Again it is enough to consider the linear stability along the principal directions. Then we 
find for g* 8 = 2g* 

f)K n i Q 9 2 M 

Sgi + 7T7 Sgi 8 - -— £ (5gi + 25g{, ), i = 1, 2, • • • , M, 



dx 5M ^ 5M 25M 2 f-r 

(223) 

# = ^^-2^g(^ + 2 ^)' -^,-,M, (224) 
r^T^, / M • l-.U • 2.---..Y. (225) 



dx 5M 

^088 _ 1 



Sgi 8 , i = M + l,M + 2,---,N. (226) 



<9x 5M 

As we can see, the fixed point is stable for any small perturbation along the principal 
directions with % = M + 1, M + 2, • • • , N. For a small perturbation along the principal 
directions of i = 1, 2, • • • , M, we can rewrite the linearized scaling equations as 

dSgi ( 3 . 2 , \ . / 2 . 4 



9a; 



("5M - 25tf ^) ^ + (oM S - - 25tf ^) ^ 



(227) 

= ( _A_ S _ A .") ^ _ _JL^ ^ <f<&. (228) 

dx V5M y 25M 2 'V ys 25M 2 lJ ys8 v ; 

Here 5y is the Kronecker delta and A,^ = 1 for all i,j = 1, 2, • • • , M. We further note that 
the combinations 5g = 5g\ + 2 <5g 88 and 5g = 2 5g 8 — 5g 88 have simple scaling equations 

()()<J ' ( ' % - ^ (229) 



<9x V5M J 25M 2 

£~55>* < 23 °> 

Note that the diagonalization of the combination g leads to one eigenvalue of — 1/M and 
M — 1 eigenvalues of 1/M. The combination g is an irrelevant perturbation, while the 
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combination g is an relevant perturbation for M > 1. Hence we conclude that only the 
one-channel S c = 3/2 fixed point is stable, while the M-channel S c = 3/2 fixed point is 
unstable for M > 1. 

Now we are going to close this section after discussing how to extend our analysis to the 
full interaction case. In this case, the exchange interaction is a little bit complicated 

H 1 = J l ® S l c S l V7 , (231) 
f = J G ® 7° + h ® «° + h ® i~° + J 68 ® 7* + J 78 <g> k* + J 8 8 ® r*. (232) 

Here Jj's are matrices whose dimension is determined by the number of the conduction 
electron partial waves. Applying the third order scaling equations, we find more fixed points. 
However, only the stable fixed points we considered in the previous section - one-, two-, and 
three-channel fixed points - remain stable in the presence of the multiple conduction electron 
partial waves. As we found in the analysis of Sec. IV, a big zoo of new unstable multi-channel 
fixed points are generated. 

Our overall conclusion is that the one-, two-, and three-channel fixed points remain stable 
with the enhanced bare exchange couplings, even in the presence of the multiple conduction 
electron partial waves. 

VI. DISCUSSION AND CONCLUSION. 

We have introduced and studied a realistic model Hamiltonian for Ce 3+ impurities with 
three configurations (/°, f 1 , f 2 ), which are embedded in cubic normal metals. Using the 
third order scaling theory, we analyzed a various terms in our model Hamiltonian. 

In our study, a new exchange interaction has been discovered, which may lead to the 
non- Fermi liquid S c = 3/2 fixed point. This model is under study using the numerical 
renormalization group (NRG) and conformal field theory |P"I|| . Unfortunately, this new 
fixed point becomes unstable against the mixing interactions between the Tg and the Tqj 
conduction electrons. 
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Within our model study, we found that the one-channel (the Tq or the Tj conduction 
electron manifold), the two-channel (the r 8 conduction electron manifold), and the three- 
channel (the r 6 ©r 8 or the r 7 ©r 8 manifold) S c — Si — 1/2 non- Fermi liquid fixed points are 
stable. From the scaling equations, we also identified several unstable fixed points. Of special 
interest are the fixed points with the higher conduction electron spins S c = 3/2,5/2,7/2, 
which may generate non-Fermi liquid physics. The multi-channel exchange interactions with 
the different conduction electron spin size are also identified at the unstable fixed points. We 
may ask what the physical meaning of these unstable fixed points are. When some specfic 
exchange interaction coupling is strong enough, the system will show the physics at high 
temperatures relevant to this exchange interaction model before flowing finally to the stable 
fixed point at zero temperature. 

The multiple conduction electron partial waves generate many unstable multi-channel 
fixed points. Again we find stability of the one-channel, two-channel and three-channel fixed 
points. 

We close by noting a slight discrepancy between our work and that of Koga and Shiba 
[TJ|, who studied the possible non- Fermi liquid ground states of an f 3 , J = 9/2, T 6 doublet 
interacting with conduction electrons. The exchange coupling with conduction electrons 
is mediated by virtual f 3 — f 2 charge fluctuations. In an NRG study of this model with 
a particular choice of coupling strengths, it was found that a fermi liquid ground state 
always arose. However, while the initially derived hamiltonian is identical in form to that 
considered in our paper for the r 6 © T 7 © r 8 sector of conduction states, we find that the 
NRG Hamiltonian used in this paper has explicit broken channel symmetry for the effective 
two-channel coupling. We believe that this may be the source of the divergent conclusions 
between our work and theirs. 
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APPENDIX A: MIXING MATRIX ELEMENTS 

In this Appendix, we shall discuss briefly how to project the atomic electron operators 
to the atomic electron Fock space. We will use the Wigner-Eckart theorem in the projection 
and calculate the fractional parentage coefficients. 

To treat the strong on-site Coulomb interaction accurately, we introduce the projected 
states belonging to the atomic electron Fock space. In our model for Ce 3+ , we assume 
the strong on-site Coulomb interaction, strong spin-orbit coupling, and strong crystalline 
electric field (CEF) in the atomic electron space. We project the atomic electron Fock 
space to the three configurations f°, f 1 , f 2 which are relevant to the trivalent Ce impurities 
embedded in normal metals. The total orbital angular momentum (L) and the total spin 
angular momentum (S) operators are good quantum numbers only in the presence of the 
Coulomb interaction. When we include the LS coupling as well as the Coulomb interaction, 
now the total angular momentum operator, J = L + S, is a good quantum number. Below 
we show the projection of the atomic electrons when the total angular momentum operator 

— * 

J commutes with the atomic Hamiltonian. The extension to the crystal electric field case is 
straightforward. 

The projection of the atomic electron operators to the atomic electron Fock space is 
realized by 

i = I/ ></°I + EI/V></VI + EI/ 2 ;j></ 2 ;^ (ai) 

U - !/#.!■ (A2) 
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Here the right hand side of 1 is the restricted completeness relation to the three configura- 
tions /°, f 1 , f 2 . The labels i,j count all the possible states in the f 1 and f 2 configurations, 
respectively. 

In the f 1 configuration, the atomic electron states are labeled by ji = 5/2, 7/2 and their 
azimuthal quantum numbers. That is, 

\f X \ni > and I/ 1 -mi > . 

In the f 2 configuration, the total angular momentum ranges from through 6. The atomic 
electron states can be represented by \f 2 LSJfi2 >• For the spin singlet states, J = L = 
0, 2, 4, 6. For the spin triplet states, L = 1, 3, 5 and J = L, L ± 1. We are now in a position 
to find the projected atomic electron operator of fj cfic . 

/icMc - l/° >< f\ fwc E l/ViPi >< /\7'iA*i| 

+ E l/^'i^i >< /.cmc E l/ 2 ^^ 2 >< z 2 ^^! 
= \f° >< f l jc»c\ + E E Hjc^fjmlfLSj^) i/Vi/ii >< ; 2 ^j/i 2 |, 

Jl/Ltl LSJfl2 

(A3) 

A0>c5 fjmlfLSJ^) = < fj^ f jotte \f 2 LSJ^2 > ■ (A4) 

Here A measures the mixing strength between J 1 and / 2 atomic states. Applying the Wigner- 
Eckart theorem, we can rewrite the mixing strength A as the Clebsch-Gordan coefficients 
multiplied by a prefactor 

AO'cMc! PjiHilfLSJfJtz) = K( Jc ; fj^fLSJ) < j c Hc,jif*i\Jf*2 > • (A5) 

Here the prefactor K(j c ; f 1 ji\f 2 LSJ) is the fractional parentage coefficient and is listed in 
Tables ^TJ,^TTT|. 

In the strong LS coupling limit, we now illustrate how we can calculate the matrix 
elements < f 1 jifi\ \ fj c ^ c \ f 2 LS : Jfi 2 > to find the fractional parentage coefficients. We 
construct symmetry eigenstates of the total angular momentum in the f 2 configuration. 
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\f 2 LS = 1 : J/i 2 > 



Y Y < 3mi;3m 2 |LM > 

M m 1 m 2 
( 



V2 



<U;LM\J^> /i lT / r L 2T |0> 
+ < 10; LM\Jfx 2 > ^ ( fl^fl 2l + fljl^ ) |0 > 
+ < 11; LM\Jfi 2 > f j mil fl 2l \0 > 



\ 



(A6) 



\f 2 LS = : J/x 2 > 



Y Y < 3m 1 ;3m 2 |LM > 



M m\m,2 



x < 00; LM\Jfi 2 > - ( /^/^ - flJU ) 1° > • 



(A7) 



These symmetry states are not energy eigenstates of the atomic Hamiltonian (Coulomb 
interaction and LS coupling) in the f 2 configuration. Mixing can occur among the same 
J multiplets, but this is neglected in the LS limit. Below we evaluate the desired matrix 
elements for the above symmetry eigenstates. 

Spin singlet case: In this case, J = L and only even L = J terms are nonvanishing. 

< fjm I /icMc I f L S = 0:Jfi 2 > 

= Y (-l) 1/2 ~ Qc < j c fi c \^a c ;3m c > < ji//i|^a c ;3mi > < Jfi 2 \3m c ; 3m! > 



•l) J+ic+3/2 V / (2j/c + l)(2ji + l) { JC 31 3 ) <J// 2 |.A-/'.-:.y l /'i>. 

3 3 1/2 



(A8) 



The final result is proportional to the correct Clebsch-Gordan coefficients. Here we in- 
troduced the 6 — j symbols . The explicit calculation gives the fractional parentage 
coefficients in terms of 6 — j symbols. The fractional parentage coefficients are listed in 
Table [XTI[ 

Spin triplet case: In this case, only odd L terms are nonvanishing. 

< fjm I fjc^c I f LS = 1 : J//2 > 

= ^ Y 12 12 < ^M|3m c ;3mi > < lm\)-a c - -on > 

ofiaic Afmmim c * * 
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x < jc^c\-jOt c \ 3m c > < jiHi\-ai, 3m 1 > < lm; LM\J^ 2 > 
= K(j c ; f 1 j 1 | f 2 LS = 1 : J) < J^jcUvh^i >, 
K{j c ;f l h | f 2 LS = 1 : J) 



6(2L + l)(2j 1 + l)(2 Jc + l) 
3 3 L 



3 jc 1/2 



x^( 2 g + i) (-i)^ 
q I j i g | | 1/2 i g 




(A9) 
(A10) 



(All) 



Here P = 2j c + 2L + 2j 2 + 1- For a given L, there are three values of J = L — 1, L, L + 1. 



The fractional parentage coefficients are listed in Table XIII. 



APPENDIX B: MULTIPLICATIVE RENORMALIZATION 

When we apply the perturbation theory to the Kondo problems, logarithmic divergences 
appear at all orders in the coupling. This infrared divergence can be handled by renor- 
malizing the exchange coupling constant for the interaction vertex, and renormalizing the 
"mass enhancement" parameter for the pseudofermion Green's function. Then we do the 
perturbation theory in the renormalized parameters in the hope that it leads to converging 
solutions. 

Though the renormalization can be defined in several different ways, we may consider 
the progressive removal of the conduction band edge states for defmiteness. At temperature 
T, only the conduction electrons (thermally excited) inside the band of order T with respect 
to the Fermi level play an important role in determining physical properties. Thus we can 
integrate out the band edge states (virtually excited states) to find the effective Hamiltonian. 
When the renormalized Hamiltonian is of the same form as the original one, we can find 
the scaling equations for the model parameters. Though the following analysis is restricted 
to the perturbative regime, we can derive qualitative results out of this. For quantitative 
results, the full numerical renormalization group method is required. 
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Here we just sketch how to generalize the multiplicative renormalization group [0,01 to 
the following form of the exchange interaction. 

Hl = J2J l ® Si(0) S}. (Bl) 

i 

Here Si = 1/2 is the impurity spin, and S c = 1/2 is the conduction electron spin. The 
exchange coupling is arbitrary square matrix. For this exchange interaction, the interaction 
vertex is 

-pa ja -pa na _ a (YK r )\ 

pqa/3;fj,u ^pr rq ^ ca.fi ^ [iv \ I 

This exchange interaction describes the most general interaction introduced in chapter II. 
Though the Schrieffer- Wolff transformation generates the potential scattering terms, we are 
not going to include them here. 



In the Kondo problem, the multiplicative renormalization is defined as |22]j23 



G Zl G, (B3) 

F -> z 2 F, (B4) 

f a -> [Z a Y l f a , (B5) 

J a -> — J a Z a . (B6) 

Z\Z2 

Here G and F are the Green's functions for the conduction electrons and the pseudo fermions 
representing the local spins. V is the scattering interaction vertex (matrix), and J is the ex- 
change coupling (matrix). Since the renormalization of the conduction electron is negligible, 
z\ — > 1. Here we note that Z a ia a matrix and a = x,y, z. 

After following the same algebraic steps as in Ref |22"1J23|| , we find the following scaling 
equations. 

Qjp- 

fd\J'my) J'my J'\rw -^f ("^inv) ■ (-^^) 

The functions fd and (matrix) are found from the irreducible or skeleton digrams of d{uS) 
and T{uj) using the perturbation theory. 
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d{u) = \ + f d {J) K(u), 

f» = l+Ff(J) K(w), 
D 



K(uj) = log 



+ iir9(w). 



(B8) 
(B9) 
(BIO) 



APPENDIX C: SCALING EQUATIONS UP TO THIRD ORDER 

We consider the following form of the interaction. 

Hi = J2J l ®S l c (0) Sj. (CI) 

i 

Here S l is the i-th component of the impurity spin operator (Sj = 1/2) and S l c (0) is the i-th 
component of the conduction electron spin operator (S c = 1/2) projected at the impurity 
site. «P is the exchange coupling matrix. 

It is straightforward to derive the scaling equations. The second order diagrams give 

-N(0)K(u) S j S l J j <g> S{J l ® SI (C2) 
N(0)K(u) S l S j J j <g> S j c f <g> S l c (C3) 

The first comes from the particle excitation, while the second from the hole excitation. Their 
contribution is 



N(0)K(u) ie ijk S k P J 1 <g> S j c Si 

= i N(0)K(u) e ijk UjijJ* ® S[ S k . 



(C4) 



The third order diagram gives 

[N{0)] 2 K{u) S k S j S l Ti[J k <g> S c fc J< <g> Sl] J j <g> S J C 

= \ [n(o)] 2 k(u) Trirr] P ® si S'S'S' 

= - [N{0)] 2 K{lo) Tt[PP] P ® S j c [ ie ijk S k S i + S j S i S i }. 
Then after the self energy correction, third order scaling equations are 



dx 

Here (i,j,k) are cyclic. 



9 i 9 k + 9 k 9 i ]- l 9 i Tr[gigt + g k g k 



(C5) 



(C6) 
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TABLES 



Table I. Crystal electric field energy eigenstates for J = 5/2 multiplet in the cubic symmetry. 



J = 5/2 Multiplet 


States 


|lf /2) ;T/i> 


- v /J|±5/2>+ v /| | T 3/2> 


|lf /2) ;+,T / j> 


| ±1/2 > 


|lf /2) ; -,T/|> 


yjl |±5/2>+yi| T 3/2> 



TABLE I. Crystal electric field energy eigenstates for J = 7/2 multiplet in the cubic symmetry. 



J = 7/2 Multiplet 


States 


ri 7/2) ;T/|> 


±75 |T7/2>± v /J|±l/2> 


|r? /2) ;t/i> 


±^ | ±5/2 > | ±3/2 > 


ri 7/2) ;+,T / j> 


± v /J| T 7/2>TV / 5l±l/2> 


\4 /2) ; -,T/|> 


±| | ±5/2 > ±^ | ±3/2 > 
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TABLE II. Decomposition of J-multiplets in cubic symmetry. All possible J-multiplets are 
listed for the f 2 configuration. 



J-multiplets Dj 


Decomposition 


D 
Di 
D 2 
D 3 
D 4 
D 5 
D 6 


r 4 
r 3 ©r 5 
r 2 © r 4 © r 5 
ri © r 3 © r 4 © r 5 
r 3 © 2r 4 © r 5 
ri © r 2 © r 3 © r 4 © 2r 5 


TABLE III. Crystal electric field energy eigenstates for J = 4 multiplet in the cubic symmetry. 


J = 4 Multiplet 


States 


irf } > 


|4> +V /J |0> +v /J|-4> 


|rf } ; + > 


^(|2>+|-2> ) 


|r 3 V> 


~4k i 4 >+75 io>- v /5i-4> 


|ri 4) ; o> 


^(|4>-|-4> ) 


|rf;±i> 


^lT3>+ V / | |±1> 


|r1 4) ;0> 


^(|2>-|-2> ) 


|rS 4) ;±i> 


v / il^ 3> -7sl ±1> 
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TABLE IV. Irreducible tensor operators for the conduction electrons. As noted in the text, 
the r 4 irrep is the spin operator. Hence there are all in all 6 possible conduction spin operators. 
The mixing between Tq and r 7 (Tq <g> r 7 = T2 © T$) does not have a spin degree of freedom. 



| r cb >a 


5 < r cb | 


Irreducible tensor operators 


r 6 s 




t 1 © r 4 


r 7 s 




ri ©r 4 


r 8 s 


$r 8 


ri © r 2 © r 3 © 2r 4 © 2r 5 


r 6 s 




r 3 © r 4 © r 5 


r 7 s 


5T 8 


r 3 © r 4 © r 5 



TABLE V. Selection rules for the Anderson hybridization between / 1 T 7 and f 2 configurations. 
For example, we can read off from this table that only the T 8 conduction electrons can mix with 
/ 1 r 7 states to reach / 2 r 3 states. 





f CEF states 


r 6 ©r 7 


r 2 ©r 5 


r 7 ©r 7 


ri ©r 4 


r 8 ©r 7 


r 3 © r 4 © r 5 
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TABLE VI. We tabulate < j c ^7', §r 7 | J2T1 >. An asterisk denotes the violation of the triangu- 
lar inequality for the addition of two angular momenta. The Ti irrep exists only in the J = 0, 4, 6 
multiplets in the f 2 configuration. 



3c 


Ji 


h 


< j c r 75 |r 7 | j 2 t 1 > 


5/2 


5/2 




4 
6 


/T 

V 3 
* 


7/2 





4 
6 


* 


TABLE VII. We tabulate < §T7| J2X2 >• The irrep is absent in the j c = 5/2 partial 
waves. The T2 irrep exists only in the J = 3, 6 multiplets in the f 2 configuration. 


jc 


Ji 


J 2 


<i c r 6; |r 7 |j 2 r 2 > 


7/2 


5/2 


3 
6 


V7/18 

1 

3 
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TABLE VIII. Table for < j c F 8 ; |r 7 | J 2 T 3 >. Note that £j 2 | < j c T 8 ; |r 7 | J 2 T 3 > | 2 = 0.5, 
which is none other than a correct normalization. The T3 irrep exists only in the J = 2, 4, 5, 6 
multiplets in the f 2 configuration. 



3c 


Ji 


h 


< j c r 8; |r 7 | j 2 r 3 > 


5/2 


5/2 


2 


ait 

V 84 






4 


FT 








V 21 






5 


1 

2 


7/2 




2 
4 
5 
6 


V 126 

V 154 

n~ 

V 18 
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TABLE IX. Table for < j c T 7 ; §r 7 | J 2 T 4 > and < j c r 8 ; §r 7 | J 2 r 4 > in Table |X|. Since the 
two r 4 irreps are present in the J = 5 multiplets, the J = 5T4 CEF states are model parameter 



dependent. However, we note that the generic hybridization form is preserved [10 



J, 



■h 



< j c r 7 ;|r 7 |j 2 r 4 > 



< j c r 8 ;|r 7 |j 2 r 4 > 



5/2 



5/2 



2 
1 
1 

75 



7/2 



1 

7 

yf2M 



12 
1 

4vTT 
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TABLE X. Table for < j c r 6j8 ; §r 7 | J2^5 >• Since the two T$ irreps are present in the J = 6 
multiplets, the CEF states are model parameter-dependent. However, we note that the generic 



hybridization form is preserved [ 10 1 



3c 


Ji 


h 


< j c r 6 ;§r 7 |j 2 r 5 > 


< j c r8;|r 7 |J 2 r5 > 


5/2 


5/2 


2 
3 
4 
5 




nr 

V 42 

nr 

V 72 

vii 


7/2 




2 


/To" 

V 81 


5 

18-/7 






3 


36 


5^ 
36 






4 


nr 

V 264 


13 
4^3^77 






5 


Vto 

18 


1 
9 






6 
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TABLE XI. We list all the possible fixed points generated from the third order scaling equations 
of our model exchange interactions. Most of them are unstable. Some of them might not exist in 
the full calculations. 



9* 


M l 


A 


Stability 


oo 


7° 


1.30 

2 2 ' 


Stable 


oo 


K ° 


1 • 3 K 

2 2^ 


Stable 


1 


r° 


1.3 J 
2 2 ' 


Stable 


2/3 


-\ (7° + t° )±^f y + | r* 


M(7° + r°) 


Stable 


2/3 






Stable 











Unstable 


1/5 


-r° i 2 t' 


3 . 5 
2 2 ' 


Unstable 


2/3 


7 ° + T° 


M(7° + t°) 


Unstable 


2/3 


K ° + T° 




Unstable 


2/11 


k° - r° ± 2 r* 


1.30,3.50 
2 2 ^ ~ 2 2 ' 


Unstable 


2/11 


7° - r° ± 2 r i 


1 . 3 ^0 , 3 . 5 

2 2 ' ^ 2 2 ' 


Unstable 


1/4 


4 ^0 _ 2 + V2 i , 4 i 
3 ' 3 3 ' ^ 3 ' 


1 • 2 ( 7° + r° ) 


Unstable 


1/4 


4 K 2 i 2V 7 ! i 4 i 
3 3 T ^ "S - K 3 T 


1 • 2 ( k° + r° ) 


Unstable 


2/35 


3 ' ^ 3 ' ^ 3 ' ^ 3 ' 


H(7° + r°) 


Unstable 


2/35 


5 , 7 i 4^ i 4 i 
3 3 3 K 3 T 


|.|( K °+r°) 


Unstable 


1 


7° + K° 


M(7° + «°) 


Unstable 


1/2 


7° + k° + r° 


1.3 (7 + K + r0) 


Unstable 


1/6 


7 ° + k° - r° ± 2t { 


I.|(7° + K°) + l-lr0 


Unstable 


1/2 


K _ 1 ( 7 + r ) ± 2^2 y + 2 T i 


i-l(7 + -° + r ) 


Unstable 


1/2 


7 _ 1 ( K + T ) ± 2^2 K i _ 2 r i 


1.3 (7 + K +r0) 


Unstable 


2/9 


K + 4 y _ 2 T ± 2^2 y + 4 Tl 


i ■ | «° + l-2 (7° + r° ) 


Unstable 


2/9 


7 + 4 K _ 2 T ± 2^2 ^ _ 4 r i 


i.3 7 o + 1.2(«o + r°) 


Unstable 


1/18 


K _ 5 y + 7 T ± 4^5 y + 4 T% 


1.3«0 + 5.7 (7 +Ta) 


Unstable 
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1/18 


7 _ 5 K + 7 T ± 4^5 j _ 4 T i 


1 3 

2 2 


7° + |-i (»° + r°) 


Unstable 


1/42 


3 7 o _ 7 K o _ 7 r o ± 2 ^ y 


7 
2 


| (7 o + K o + r0) 


Unstable 












1/42 


| 7 o + 3 K o 7 T o ± 2^5 y 


7 
2 


|( 7 o + K o + T o ) 


Unstable 




±2^/3 k* - g t* 









TABLE XII. Fractional parentage coefficients for K(j c ;f 1 ^ | f 2 LS = : J). * means the 
triangular sum is violated. Note that for a given (S, L, J), sum of the squared fractional parentage 
coefficients over all possible j c = 5/2,7/2 is always 6/7. 



(S,L,J) 


3c = 5/2 


jc = 7/2 


(0,0,0) 


/ef 

V 7 


* 


(0,2,2) 


6 
7 


nr 

V 49 


(0,4,4) 


/22~ 
V 49 


/20 
V 49 


(0,6,6) 


* 


/ef 

V 7 
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TABLE XIII. Fractional parentage coefficients for K{j c \f\ \ f 2 LS = 1 : J). * means the 
triangular sum is violated. Hund's rule ground multiplets correspond to (S = 1, L = 5, J = 4). 
Note that for a given (S,L), sum of the squared fractional parentage coefficients over all possible 
j c = 5/2, 7/2 and over all J = L, L ± 1 is always 20/7. 



(S,L,J) 


jc = 5/2 


ic = 7/2 


(1,1,0) 


/~8 

V 7 


* 


(1,1,1) 





1 


(1,1,2) 


V 49 


V 49 


(1,3,2) 


/54 
V 49 


/l6 
V 49 


(1,3,3) 





1 


(1,3,4) 


nr 

V 147 


V 147 


(1,5,4) 


/220 
V 147 


V 147 


(1,5,5) 





1 


(1,5,6) 


* 


IT 

V 7 
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FIGURES 

FIG. 1. Schematic display of all the possible exchange interactions. With the impu- 
rity pseudo-spin 5/ = 1/2 for /V = 5/2r 7 doublet, there are in total 6 exchange interactions 
for r6,r7, and Tg conduction electrons. In addition to the well-known one-, and two-channel 
S c = 1/2 exchange interactions with antiferromagnetic (AFM) coupling, there arise one-channel 
S c = 3/2 exchange interactions with antiferromagnetic coupling (not known in the literature), and 
one-channel S c = 1/2 exchange interactions with ferromagnetic (FM) coupling. There are also 
the mixing exchange interactions through which the conduction electrons can change their CEF 
symmetry labels after interacting with the impurity pseudo-spin. 

FIG. 2. The scaling diagrams up to the third order. The solid lines denote the conduction 
electron propagator and the dotted lines the impurity spin. The external solid lines are confined 
to the conduction electrons close to the Fermi level, while the intermediate solid lines denote the 
conduction electrons virtually excited to the high energy states. At each vertex, the exchange 
coupling matrices are determined by the specific exchange interaction type. When the multiple 
conduction electron partial waves are considered, the exchange coupling matrices should be written 
as the direct product of the channel matrix and the exchange interaction matrix. 

FIG. 3. Schematic screening of Si = 1/2 by S c = 3/2. In the first shell screening, the 
energy is minimized by the compensation of the impurity spin by the conduction electrons of 
S c = —3/2, —1/2. This picture does not violate the Pauli exclusion principle. After the first shell 
screening, the effective impurity spin is Si = 3/2. In the second shell, S c = —3/2, —1/2 conduction 
electrons cannot come closer to the impurity site due to the Pauli exclusion principle. Hence the 
effective exchange coupling becomes antiferromagnetic. In this case, the effective impurity spin is 
Si = 1/2. Thus the effective impurity spins alternate between Si = 1/2 and Si = 3/2 with onion 
shell screening process. This alternating picture will be verified in Chapter V using the numerical 
renormalization group technique. 
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FIG. 4. Flow diagram in the Tg conduction electron space. There are three stable fixed 
points in the Tg conduction electron space. Only the two-channel fixed point remains stable when 
the full exchange interactions are included. The other two fixed points derived from S c = 3/2 
conduction electrons become unstable against the mixing exchange interactions between the Tqj 
and the Tg conduction electrons. See the text for more details. 
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